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ABSTRACT 

We describe two related experiments that measured the anisotropy in the cosmic microwave background 
(CMB). QMAP was a balloon-borne telescope that flew twice in 1996, collecting data on degree angular 
scales with an array of six high electron mobility transistor-based amplifiers (HEMTs). QMAP was the 
first experiment to use an interlocking scan strategy to directly produce high signal-to- noise CMB maps. 
The QMAP gondola was then refit for ground based work as the MAT/TOCO experiment. Observations 
were made from 5200 m on Cerro Toco in Northern Chile in 1997 and 1998 using time-domain beam 
synthesis. MAT/TOCO was the first experiment to see both the rise and fall of the CMB angular 
spectrum, thereby localizing the position of the first peak to Ipeak — 216 ± 14. In addition to describing 
the instruments, we discuss the data selection methods, checks for systematic errors, and we compare 
the MAT/TOCO results to those from recent experiments. We also correct the data to account for an 
updated calibration and a small contribution from foreground emission. We find the amplitude of the 
first peak for 160 < I < 240 to be STpeak = 80.9 ± 3.4 ± 5.1 /iK, where the first error is statistical and 
the second is from calibration. 



1. INTRODUCTION 

Experiments aimed at measuring the anisotropy in 
the CMB require a combination of sensitive detectors 
and novel observing strategies. The observational goal 
is to measure micro-Kelvin celestial variations in ther- 
mal emission with a telescope observing from an envi- 
ronment that is some ten million times brighter. Be- 
low 90 GHz, the detectors of choice have been high 
electron mobility transistor based amplifiers designed at 
NRAO [HEMTs, Pospieszalski (1992)]. Above 90 GHz, 
bolometers are the best detectors [e.g., (Bock et al. 1998; 
Lee et al. 1996; Downey et al. 1984; Tucker 1991)]. SIS- 
based systems (Kerr et al. 1993) near 100 GHz have 
the speed and intrinsic sensitivity of transistor am- 
plifiers though do not yet have the large instanta- 
neous bandwidth of bolometers or HEMTs. Over the 
past five years, instruments have been designed for di- 
rect mapping of the CMB [e.g., QMAP, BOOMERanG 



(Piacentini et al. 2001), MAXIMA (Hanany et al. 2000), 
TopHat (2001) ], and for beam synthesis [Saskatoon 
(SK) (WoUack et al. 1997), MSAM (Fixsen et al. 1996), 
Tenerife/Bolo (Romeo et al. 2000), PYTHON 
(Coble et al. 1999), VIPER (Peterson et al. 2000)]. 
More recently, interferometers based on HEMT am- 
plifiers have reported CMB anisotropy results [CAT 
(Baker et al. 1999), DASI (Leitch et al. 2001), lAC 
(Harrison et al. 2000), CBI (Padin et al. 2000)]. Though 
the primary data product of the interferometers and beam- 
synthesis experiments is the angular spectrum, data taken 
with these techniques can be turned into maps (e.g., 
Tegmark 1997). Conversely, there is always some filtering 
involved in the mapping experiments. A common element 
of these experiments is that they are limited by systematic 
error. 

In this paper, we describe the instruments for the 
QMAP experiment and the Mobile Anisotropy Telescope 
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on Cerro Toco (MAT/TOCO or TOCO for short). We 
supply the details necessary for assessing the quality of the 
data and reproducing the experimental method. QMAP 
is described in part in Devlin et al. (1998) and analyses of 
the data are presented in Herbig et al. (1998), de Oliveira- 
Costa et al. (1998, 1999), Xu et al. (2000), and Park et 
al. (2001). Balloon borne mapping experiments have a 
long history (Weiss 1980; Partridge 1995) though highly 
interlocking scan strategies over limited regions of sky are 
more recent (Staren et al. 1999; de Bcrnardis et al. 2000; 
Lee et al. 2001). QMAP, which flew twice in 1996, was 
the first of these to produce a "true map" of the CMB, 
complete with pre- whitening and full covariance matrices. 
QMAP was comprised of a focal plane array of three dual 
polarized HEMT channels with an angular resolution of 
roughly 0?8. The beam array was steered on the sky by a 
large chopping flat. 

TOCO used the QMAP gondola and receiver refit to 
work with a mechanical cooler instead of liquid cryo- 
gens. It also employed two SlS-based^^ 144 GHz detec- 
tor systems to improve the resolution to 0?2. TOCO 
employed the Saskatoon-style beam synthesis strategy 
(Netterfield et al. 1995) with eight independent detec- 
tors. Instead of observing near the NCP from Saska- 
toon, Canada, we observed near the SCP from the side 
of Cerro Toco in Northern Chile At 144 GHz, the 
atmospheric column density in Saskatoon is too large 
for anisotropy measurc;nients; a high altitude site such 
as the Chilean Altiplano is required. TOCO operated 
for two seasons in 1997 and 1998. The primary re- 
sults and short description of the instrument are given 
in Torbet et al. (1999) and Miller et al. (1999). The 
0?2 resolution allowed us to locate the first peak in the 
angular spectrum at I « 212 ± 14 (Knox and Page 
2000). In the context of the popular adiabatic CDM 
models, this shows that the universe is geometrically 
flat (Doroshkevich et al. 1978; Kamionkowski et al. 1994; 
Bond et al. 1994; Hu & White 1996; Cornish 2000). 

2. OVERVIEW OF GONDOLA AND MOBILE ANISOTROPY 
TELESCOPE 

The TOCO experiment is shown schematically in Fig- 
ure 1^^. The part of the flgure containing the optics and 
receiver is the QMAP balloon gondola. The radar trailer, 
on which the gondola is mounted, has a separable mag- 
nesium base. Three legs hold the base off the trailer and 
stabilize it. When the telescope is transported, the legs are 
removed. The radar trailer has a 1.4 m diameter precision 
bearing on which is mounted a 2.5 cm thick flat plate that 
holds the gondola. The plate and gondola are rotated us- 
ing an on-axis DC motor^*. CMB observations are made 
with the telescope in a fixed position; a brake holds the 
telescope in place allowing the motor to be shut off during 
observations and preventing it from drawing large currents 
or oscillating as it seeks the target position in high winds. 
The motor has a 15 cm diameter hole in the center through 
which cables and refrigerator hoses pass from the inside of 

SIS stands for Supcrconductor-Insulator-Supcrcoiiductor. The detecting clement is a quasi-particle mixer (Tucker & Feldman 1985). 
The Cerro Toco site of the Universidad Catolica de Chile was made available through the generosity of Prof. Hernan Quintana, Dept. of 
Astronomy and Astrophysics. It is near the ALMA site. 

1^ The Nike Ajax radar trailer was donated to the University of Pennsylvania by Lucent Technologies. 
1® The motor is a Compumotor DR 1100 A - 100 Nm torque. 



the telescope to the outside. The compressor that runs 
the mechanical cryocooler is mounted on the trailer. The 
azimuth is instrumented with an absolute 17-bit encoder 
and a 20-bit resolver. 

3. THE RECEIVER 

Radiation from the sky enters the dewar through a 15.25 
cm diameter vacuum window made of 0.56 mm polypropy- 
lene and is collected with corrugated feed horns as shown 
in Figure 2. Three aluminum baffles deflne the entrance 
aperture, one is attached to the 40 K cold plate, one is at- 
tached to the dewar just inside the vacuum window, and 
one is attached to outside of the dewar. Strips of alu- 
minized Mylar connecting the cold feeds to the ambient 
temperature dewar block RF interference and reduce op- 
tical loading on the cold stage. To prevent the formation 
of frost on the window, warm air is passed through a vol- 
ume in front of the vacuum window defined by a Saran 
Wrap-covered aluminum cone attached to the outside of 
the dewar. 

Two NRAO SIS mixers are attached to the w 4.5 K 
stage of the dewar and six HEMTs are attached to the 40 K 
stage, two with center frequencies of 31 GHz (in Ka band) 
and four with center frequencies of 42 GHz (in Q band). 
Warm amplifiers, bandpass filters, noise sources, and local 
oscillator for the SIS system are housed in a 293 K "back- 
pack" attached to the outside of the dewar. The primary 
difference between the TOCO and QMAP receiver config- 
urations is that QMAP used liquid cryogens and TOCO 
used a mechanical refrigerator to cool the HEMTs as well 
as the SIS mixers. 

3.1. The SIS Configuration 

The SIS configuration is shown in Figure 3. The mix- 
ers are coupled to the sky with conical corrugated feed 
horns. A round-to-rectangular transition at the base of 
the horn transforms the incoming signal to a single polar- 
ization which is added to a signal from the local oscillator 
using a 20 dB branch line coupler. The combined signals 
are fed into the SIS mixer block. A sliding backshort, con- 
nected to the outside of the dewar by a flexible shaft, is 
used to tune the SIS mixer, while cold, for optimal cou- 
pling. The SIS is biased through a bias T which allows 
transmission of the RF signal while blocking the DC bias 
voltage. An associated circuit controls and reads the cur- 
rent through and voltage across the SIS. The signal then 
comes out of the bias T, goes through a 20 dB 4-6 GHz 
isolator (P & H Labs) and into a cold C-band (3.95-5.85 
GHz) HEMT amplifler with 33 dB gain, as diagrammed 
in Figure 3. The SIS, branch line coupler, bias T, and 
C-band HEMT amplifier were all designed and fabricated 
at the National Radio Astronomy Observatory (NRAO). 

The output of the cold HEMT is carried on stainless 
steel semi-rigid coax from the cold stage to the backpack 
where it is amplified by a warm 44 dB RF amplifier, fil- 
tered through a 3.5-6.5 GHz filter and detected with a 
detector diode (Hewlett-Packard: ±0.2 dB flatness and a 
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Fig. 1. — Schematic of the Mobile Anisotropy Telescope. The legs are shown in their stowed position on the left end of the trailer. When 
observing, the compressor is also on the trailer. The QMAP gondola comprises everything above the bearing. 




Fig. 2. — View of the inside of the dewar. Shown are the 4.5 K mounting structure, the thermal shield and cold plate at 40 K, and the 
shield at Ri 200 K. In the QMAP configuration, the entire cold plate was cooled in the lab to 4.2 K with liquid Helium, or to 2.7 K in flight 
(~ 33, 000 m altitude). Pressure regulated LN2 held the intermediate temperature stage near 77 K. The SIS mixers are mounted to the 4.5 K 
structure. There is a small section removed from the bottom of the chopper to accommodate the outer feed baffle. One source of modulated 
radiation is the cavity formed by the moving chopper and the feed baffle. 




Fig. 3. — Components of the SIS Receiver. The components on the right are in the dewar. Components on the left are housed in the 
thermally stabilized 293 K backpack. The temperatures of all critical components are monitored. 



typical sensitivity of 300 mV/mW). To increase the lin- 
earity and to bring the detector output into a convenient 
range (« 0.01 V), the diode is shunted with a ~ 1 kfi 
resistor. The RF filter ensures that only RF radiation in 
the passband of the C-band IF amplifier makes it through 
the system. The detected output is then amplified by 100 
and sent to the data acquisition system. The IF system 
temperature was measured in the lab with a beatable 50 il 
load and found to be TgyslIF) =< 6K. The net gain of the 
system is 73 dB. 

The LO (Millimeter Wave Oscillator Company) consists 
of a cavity stabilized oscillator with an output frequency 
of 48 GHz which injection-locks a second oscillator, which 
in turn injection locks a high power Gunn oscillator. The 
oscillator feeds a tripler with an output frequency of 144 
GHz. This signal then goes through an attenuator, into a 
3 dB directional coupler and is carried by WR-6 waveguide 
to the cold stage. On the cold stage a 90° hybrid sphts 
the LO power and sends it to the two mixers. The third 
port of the coupler is used to couple test and calibration 
signals into the system. The total LO power is about 8.5 
mW, ~ 300 nW of which impinges on each SIS mixer. 

3.2. The HEMT Configuration 

Figure 4 shows the configuration of the HEMT re- 
ceivers for the TOCO 1997 and 1998 observing seasons. 



The QMAP configuration is almost identical. The ra- 
diation enters the dewar through the vacuum aperture 
described above and is collected using Ka- and Q-band 
corrugated horns. An orthomode transducer (OMT) 
splits the radiation into vertical and horizontal linear po- 
larizations. The Ka-band OMT has a reflection coef- 
ficient of less than —26 dB over the bandpass of the 
channel. The Q-band OMTs have a reflection coeffi- 
cient of less than —20 dB. Each polarization is carried 
by rectangular copper waveguide to a separate low noise 
NRAO cryogenic HEMT ampUfier(Pospieszalski 1992; 
Pospieszalski et al. 1994). The passband of the system is 
set by the amplifiers and waveguide cutoff at the low end 
and a filter at the high end. Bends in the waveguide are 
limited to a radius of curvature of at least 1.5 cm in order 
to minimize reflections (Monnelly 1996). 

At the input to each HEMT amplifier is a 20 dB crossed- 
guide Bethe hole coupler that is used to inject a noise 
pulse. There is one source for all four of the Q-band chan- 
nels and one for both Ka-band channels. 

Upon leaving the dewar, each of the Ka-band signals 
travels through waveguide to a warm amplifier. Ka2, the 
vertically polarized channel (Table 1), is amplified with a 
single warm amplifier with a gain of about 50 dB and a 
noise figure of 3-4 dB {Tsys « 400 K). The signal from Kal, 
the horizontally polarized channel, was found to be com- 
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Fig. 4. — Components of a Q-band receiver, one of six HEMT radiometer channels. Components on the right are mounted to the 
40 K (for TOCO, 2.7 K for QMAP) plate inside the dewar, and components on the left are in the 293 K backpack. The HEMT 
output is carried on a rectangular copper waveguide to the edge of the cold stage where it is connected to 0.5 m long thin walled 
stainless steel waveguide that runs down the vacuum chamber inside the dewar to warm stage. The stainless steel waveguide 
provides thermal isolation between the cold stage and ambient temperature. Copper waveguide from the ambient temperature 
end of the stainless steel segment runs to the vacuum seal which is formed with a piece of 0.013 mm thick kapton tape. 



pressed in a similar configuration so it is passed through 
an amphfier (+29 dB), then through an attenuator (-12 
dB) and through a second amplifier (also with +29 dB) 
before being detected with a diode. The net gain of the 
Ka receiver chains is ~ 72 dB. 

The Q-band signals, upon leaving the dewar, are carried 
via copper rectangular waveguide to 1 dB insertion loss fil- 
ters (Spacek) that filters out the low frequency tail of the 
60 GHz atmospheric oxygen line. Following the filters the 
signal is amplified, sent through a band leveler (designed 
individually for each channel by Pacific Millimeter) and 
amplified again. The amplifiers are connected to the band 
leveler with a K-connector to waveguide transition. Each 
amplifier has a gain of 26-28 dB and a noise figure of « 4.5 
dB. Tests of the system with and without the band level- 
ers are discussed in Monnelly (1996). The net gain of each 
Q-band system is ~ 75 dB. 

A diode detector (Millitech) at the output of the final 
amplification stage of the HEMT channels converts the 
incident power to a voltage. The diode sensitivity varies 
over the passband and with temperature but is typically 
w 2000 mV/mW. The output of the diode is connected to 
a low noise pre-amp with a typical gain of 100 and roughly 
2 kfl input impedance. We tune the DC output by select- 
ing this impedance. Values vary by about a factor of 100 
between different channels. The output of the pre-amp is 
buffered and sent to the data acquisition system. 

3.3. Thermal, Mechanical, and Magnetic Considerations 

In the QMAP configuration, the HEMT amplifiers and 
the SIS mixers are heat sunk to a cold plate that forms 
the bottom of the liquid helium reservoir. The outer tank 
of the dewar is filled with liquid nitrogen. The vapor pres- 
sures of both the nitrogen and helium are held at a con- 
stant value with mechanical regulators. The combination 
of pressure regulated cryogens and a balloon borne dewar 
produces a constant-temperature cold plate with minimal 
microphonics. 

Due to the expense and difficulty of transporting liquid 
cryogens to Cerro Toco, the TOCO receiver is cooled with 
a Sumitomo cryogenic refrigerator (SRDK-408BA). The 
40 K stage has 40 W of cooling power. The loading on 
the 40 K stage is dominated by HEMT LED's, waveguide 
connections between the cold and warm stages, and opti- 
cal loading through the window. From the liquid cryogen 
boil-off rate in the QMAP configuration, we determined 
that the total loading on this stage is « 20 W. The total 
loading on the 4 K stage, where there is 1 W of cooling 
power, is « 400 mW. 

In Chile, diurnal variations in the optical loading led 



to variations in the temperature of the SIS mixers of 
« 400 mK. However, for the purposes of CMB anisotropy 
analysis, we used only data taken during the night when 
the temperatures were stable to better than 50 mK. The 
temperatures were monitored continuously and the small 
drifts accounted for in the analysis through a calibration 
model. There was also a temperature variation of approxi- 
mately 10 mK synchronous with the refrigerator drive mo- 
tor operating at 1.2 Hz. This term is asynchronous with 
the chopping frequency, and the resulting gain fluctuations 
are too small to affect our results (Section 13). 

The backpack that houses the warm electronics mounts 
on the dewar and is thermally controlled to within ±0.5 K 
over a typical night of observations. As with the cold stage, 
the temperature is monitored continuously and the effects 
of drifts on gain are corrected in software. All microwave 
amplifiers are heat sunk to a common aluminum plate and 
their temperatures are monitored. 

The enclosure for the electronics has two levels of RF 
shielding and is filled with pieces of microwave absorber 
wrapped in plastic bags that serve the dual purpose of 
thermal insulator and absorber of stray microwave radia- 
tion (perhaps from imperfect joints). All waveguide joints 
are wrapped with an absorber to prevent leakage of ra- 
diation into the rest of the system. We did not find any 
evidence for correlations between channels due to the in- 
strument. 

If changes in RF impedance due to vibrations are syn- 
chronously modulated by the motion of the chopping mir- 
ror, the resulting microphonic lines can mimic a celestial 
signal. Such "microphonics" can couple into the data, for 
example, through motion of the feeds or through strain in 
the microwave joints. If the coupling is large and variable, 
it can also affect the data even if it is asynchronous with 
the chopper. 

In order to minimize this coupling, the cold head motor 
is vibrationally isolated from the dewar (which is bolted 
to the gondola frame) with a set of compensating flexi- 
ble vacuum bellows. We also use bags of #9 lead shot 
placed on the cold head to damp vibrations. The con- 
nections between the cold head and the cooled electronics 
are made with strips of high purity flexible copper braid 
that efficiently conducts heat and vibrationally isolates the 
detectors from the refrigerator head. 

Although microphonic levels were low at the beginning 
of the 1997 season, a microphonic coupling developed over 
the campaign that rendered the D-band data unusable for 
CMB observations. The coupling was manifest as 1.2 Hz 
(the cold head cycle) wings of a broad 90 Hz line suggest- 
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ing amplitude modulation of a 90 Hz vibrational line. The 
souree of the vibration was traced to a combination of the 
azimuthal drive motor and the chopper. It was corrected 
for the 1998 season by modifying these two systems. The 
electronic interference from the azimuthal drive motor was 
eliminated by installing a brake. The brake allows the mo- 
tor to be shut off during CMB observations and prevents 
large currents from being drawn by the motor working to 
counteract wind loading. In addition, the chopper-induced 
vibrations were reduced by replacing the bearings in tlic^ 
chopper with flex pivots. As a result, microphonic levels in 
the TOC098 D-band data were negligible. Microphonics 
were not a problem either year in the HEMT data. 

The Josephson junctions in the SIS mixers are sensitive 
to magnetic fields. Helniholz coils aromid the dewar were 
used to measure the dependence. With the coils absent, 
the area was mapped with a Gauss meter to ensure that 
the AC fields from the chopper drive and cold head motor 
would not contaminate the data. To minimize potential 
magnetic coupling, high magnetic permeability material 
(mu-metal) was wrapped around the outside of the chop- 
per coils and around the cold head. Not only is the mag- 
netic field negligible, but the chopper synchronous compo- 
nent corresponds to values of / that do not enter into the 
CMB analysis. 

3.4. Receiver Characteristics 

The QMAP/TOCO receiver was characterized in the 
lab before each campaign, but the most rek;vant charac- 
terizations are done while observing. We use the following 
definitions: 



/ vg{v)dv 
j g{v)dv 



S g{v)di 



(1) 



The SIS is operated in double side band (DSB) mode. 
By convention we use the intermediate frequency (IF) 
noise bandwidth in the radiometer equation and report 
"double sideband noise temperatures," Tdsb, because our 
source fills both RF bands. If a source fills just one RF 
band, a "single sideband noise," Tssb, is reported. In 
an ideal system, Tdsb = Tssb/^- We use the mean of 
the USB and LSB bandwidths for the noise bandwidth in 
the radiometer equation. A full calculation of the noise 
includes contributions from the mixer and the IF ampli- 
fier (Blundell 1992; Kerr et al. 1997); for our purposes we 
treat these as lumped elements. 

3.4.1. SIS Sensitivity 

Measurements of the SIS sensitivity have been made in 
several configurations as shown in Table 2. The values of 
Tree from laboratory measurements are better than those 
made in the field. To investigate this discrepancy, we built 
an external cold load that could be bolted onto the front 
of the receiver to mimic the sky. We found that it was pos- 
sible to approximately reproduce the system temperatures 
measured with the internal load using an external load pro- 
vided that the tuning parameters were readjusted. This 
retuning compensates primarily for the change in temper- 
ature of the mixer. In the 1998 season, the combination 
of thermal loading and increasingly poor refrigerator per- 
formance led to an increase in the SIS temperatures and a 
temperature distribution different to that in the lab. Even 
though the system was tuned in the field, the lab perfor- 
mance was not duplicated. 

The SIS system noise exhibits a l/f characteristic pre- 
sumably due to the C-band HEMT, though this has not 
been verified. The l/f noise is parametrized following 
equation 2 in the last two lines of Table 2. 



where g is the receiver passband, Vc is the effective center 
frequency and A„i/ is the noise bandwidth (Dicke 1946) 
for the radiometer equation. 

The noise of HEMT -based amplifiers has a 1// charac- 
teristic (Jarosik 1996; WoUack 1995). To parametrize it, 
we fit the power spectrum of the detector diode output to 
the following form: 



<-sys 



1 



tAuV 
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3.4.2. HEMT Sensitivity 

The same set of six HEMT amplifiers were used for the 
two QMAP flights and the two observing seasons of the 
TOCO experiment; four of this set were used for the Saska- 
toon measurement. They were tested in the lab before each 
set of observations and the sensitivities were analyzed for 
each data set (Table 3). There is evidence of degrada- 
tion in the HEMT performance between the QMAP and 
TOCO campaigns above that which is expected due to the 
difference in body temperature (Pospieszalski 1989). Gen- 
erally HEMT system noise is expected to increase roughly 
one Kelvin for each Kelvin of increased ambient temper- 
ature. We suspect slow deterioration in the mechanics of 
the microwave/bias structures over the hundreds of cy- 
cles and sometimes rough handling these amplifiers expe- 
rienced. The chips were unpassivated InP so there may 
some deterioration in the chip performance though this 
has not been confirmed. 

4. OPTICS 

The telescope optics are similar to those used in the 
Saskatoon experiment (Wollack et al. 1997). Corrugated 
feeds underilluminate a 0.86 m primary mirror which in 
turn under illuminates a flat (1.8 mxl.2 m) chopping mir- 
ror (chopper). Each of the eight channels detects a sin- 
gle mode of a diffraction limited beam. The chopper is 

Throughout the paper, we use for RF frequencies and / for audio (< 20 kHz) frequencies. 



(2) 

where T is the system sensitivity in units of K-sec^/^, t 
is the integration time, and AG/G is the fractional gain 
fluctuation that gives rise to the l/f form. When the 
bandwidth is very large, the gain fluctuations dominate 
as shown by Wollack & Pospieszalski (1998). Note that 
a 1// noise spectrum corresponds to a = —1 in the vari- 
ance ("power") of the detector output. We define the l/f 
knee, fknee^ to he the frequency where the power spectrum 
(square of equation 2) increases by a factor of two over the 
value at high frequencies. 

The instrument bandpasses, giv), are measured in the 
lab for each channel. Table 1 lists the center frequencies 
and effective noise bandwidths for each channel for each 
observing campaign as calculated from equation 1. 
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Table 1 

Center frequency and noise bandwidths for each campaign. 



QMAP TOCO 97 TOCO 98 





V, (GHz) 


Ai/ (GHz) 


V, (GHz) 


Ai/ (GHz) 


V, (GHz) 


Ai/ 


Kal 


31.7 ±0.3 


4.7 ±2.2 


31.7 ± 0.3 


4.7 ± 2.2 




4.7*'^ 


Ka2 


30.8 ± 0.2 


6.2 ±0.3 


30.8 ±0.2 


6.2 ± 0.3 


32.0 ±0.1 


8.8 ±0.1 


Ql 


41.4 ±0.2 


6.9 ±0.3 


41.1 ±0.2 


6.9 ±0.3 




6.9* 


Q2 


41.3 ±0.2 


7.5 ± 0.3 


41.0 ±0.2 


7.3 ± 0.3 




7.3* 


Q3 


42.1 ±0.2 


6.3 ±0.3 


41.7 ±0.2 


6.3 ±0.3 


41.5 ±0.1 


5.3 ±0.1 


Q4 


41.2 ±0.2 


7.0 ± 0.3 


41.1 ± 0.2 


4.6 ±0.3 


41.2 ±0.1 


4.5 ±0.1 


Dl (USB) 






149.2 ±0.2 


2.9 ±0.2 


148.7 ±0.1 


2.3 ±0.1 


Dl (LSB) 






138.3 ±0.2 


3.2 ±0.2 


139.4 ±0.1 


2.6 ±0.1 


Dl (DSB)'' 






141.8 ±0.7 


5.6 ± 0.4 (3.1) 


143.7 ±0.1 


4.9 ±0.1 


D2 (USB) 






148.7 ± 0.7 


2.5 ± 0.4 


148.6 ± 0.1 


2.4 ±0.1 


D2 (LSB) 






138.9 ± 0.7 


1.7 ±0.4 


139.4 ±0.1 


2.8 ±0.1 


D2 (DSB) 






145.4 ± 3.0 


3.6 ±1.3 (2.1) 


143.3 ± 0.4 


5.2 ±0.1 



"The * indicates that the bandwidth is assumed from the previous year. 

''The D-band channels were not used in QMAP so no value is given. The SIS bandpasses were re-measured between 1997 and 1998 because they 
change with SIS tuning. The numbers in parentheses following the full RF bandwidths is the effective IF bandwidth for noise calculations. 



Table 2 
sis system parameters 



Dl (lab int.) Dl (field) (lab int.) (field) 

Tphys(K)'' 4.3 4.9 4.4 4.9 

VsimMf 13.4 13.8 14.2 14.6 

fosBivoK^fsf 0.6 1.3 0.7 1.2 

Tsy/ 30 65 ••• 61 

Tree" 26 48 ••• 44 

Trec(y-factor)f 27 ••• 35 

a ■■■ -1.0 ••• -0.8 

/fcnee (Hz) • • • 18 • • • 12 



"The physical temperature of the SIS body. The uncertainty is « ±0.2 K. At T > 5 K the SIS sensitivity is markedly reduced. 
''The optimal SIS bias voltage, Vb, across the 6 SIS junctions. 

°The total power DSB sensitivity computed from the noise power spectrum at 200 Hz, where atmospheric fluctuations are negligible, and 
the responsivity. These are for a Rayleigh-Jeans source. In the field, the atmosphere and telescope contribute f« 17 K. In the lab, the load 
contributes Ri 4 K. The loss from the feeds is measured to be negligible. 

"^Thc Raylcigh-Jcans system temperature computed from measured sensitivity and the noise bandwidth. 
''The Rayleigh-Jeans receiver temperature computed from T^ys- 

■'^The Rayleigh-Jeans receiver temperature measured with a variable temperature cryogenic load. 



a resonant, computer controlled mirror, that scans in the 
azimuthal direction while the rest of the optics remain 
fixed in azimuth and elevation as the sky rotates through 
the beams. The telescope sits inside an aluminum ground 
screen which is fixed with respect to the primary mirror, 
the receiver, and the chopper mount. 

4.1. The focal plane 

Conical corrugated feed horns receive radiation from the 
sky and transform the incident fields so they may prop- 
agate through waveguide. All of our feeds were fabri- 
cated by Custom Microwave from electroformed copper 

20 In Ka band, the depths of grooves 1-10 are: 0.411, 0.368, 0.351, 0.335 
by Ed WoUack. 



over an aluminum mandrel. They arc gold coated to sta- 
bilize the surface. The general electromagnetic design fol- 
lows the guidelines in Clarricoats and Olver (1984) and 
MacA Thomas (1978). The throat section, where the cor- 
rugations adiabatically transform from A/2 depth to A/4 
depth as the hybrid mode detaches from the feed wall, was 
designed by Wollack (1994) following the work of James 
(1982) 20 The VSWR is less than 1.05 across the waveguide 
band and the loss in the feed is negligible. We model the 
feeds with a commercial code (YRS Assoc.), CCHRN, that 
solves for the full electromagnetic field that propagates in 
the feed. Table 4 summarizes parameters of the feed horns 

, 0.323, 0.312, 0.302, 0.292, 0.282, and 0.274 cm. This feed was designed 
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Table 3 

hemt amplifier system parameters 



QMAP^- TOCO 97*^ ^ TOCO 98 





T 


T 


a 


fknce 




T 




T 


a 


fknce 




K 


luK 




IIz 


K 


iiiKv's 


K 


luK v's 




Hz 


Kal 


22 


(0.40) 0.36 


-1.2 


62 


89 


1.3 


162 


(2.6) 2.4 


-0.90 


14 


Ka2 


23 


(0.42) 0.32 


-0.92 


141 


63 


0.8 


59 


(0.9) 0.62 


-0.85 


95 


Ql 


17 


0.21 






91 


1.1 


84 


(1.2) 1.0 


-0.84 


64 


Q2 


22 


(0.60) 0.25 


-0.82 


686 


145 


1.7 


114 


(3.3) 1.3 


-1.57 


168 


Q3 


155 


(2.1 ) 2.0 


-0.78 


20 


63 


0.8 


80 


(0.9) 0.72 


-0.84 


37 


Q4 


53 


(0.85) 0.64 


-0.71 


67 


156 


2.3 


88 


(1.6) 1.3 


-0.88 


42 



"From the first QMAP flight as shown in (Dcvhn et al. 1998). Ql did not work during the flight so wc report the lab measurements 
(Monnefly 1996). The four entries for each campaign correspond to the system temperature, the measured system sensitivity at 100 Hz 
(in parentheses) and very high frequencies, the gain fluctuation exponent, and the 1// knee. In the fits to the power spectra, / < 5 Hz is not 
included. Because of the 1// HEMT characteristics, T = Tsys / \/ is always smaller than the measured value of T at 100 Hz. With the 
centroids in Table 1, the full noise spectrum may be recovered. 

''Prom measurements in the field at 100 Hz. 



Table 4 

Design Parameters for All Feeds 





Ka 


Ql/2 


Q3/4 


D 


Scmiflarc angle 6o (dcg) 

Skyward aperture diameter d^(cm) 

OMT aperture diameter (cm) 

Beamwidth 9lZT (deg) 

Phase error A (MacA Thomas 1978) 

VSWR 

Forward gain (dBi) 
Number of corrugations 
Length (cm) 


6 

4.2 

0.833 

18 

0.11 

1.03 

20.9 

71 

19.3 


4.4 
2.0 
0.650 
18 

0.065 

1.03 

20.6 

67 
12.0 


5.4 
2.1 
0.650 
16 

0.051 

1.03 

21.4 

85 
14.5 


5.4 
0.89 
0.173 
17 

0.011 

1.05 

20.8 

82 

3.96 
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used in the QMAP and TOCO experiments. 

The Ka and Ql/2 feeds were also used in the Saska- 
toon experiments where they had three concentric choke 
grooves on the flange that defined the skyward aperture. 
It was bcheved that these would reduce "edge currents" 
resulting in a cleaner beam. To move the feeds as close 
together as possible in the focal plane, we removed those 
grooves and saw no degradation in performance. 

Figure 5 shows a map of the beam pattern made by ob- 
serving Jupiter. The feed horns are arranged so that the 
D-band beams with 0.2° resolution are placed as close to 
the center of the focal plane as possible. Beam parame- 
ters for each channel have been calculated and measured 
for each campaign as shown in Table 5. 

Uncertainty in the beam solid angle leads to calibra- 
tion uncertainty. Confidence in our understanding of the 
beams has been built by observing a variety of sources over 
two years and comparing those results with detailed mod- 
els. QMAP mapped the beams using Cas-A and Saturn; 
TOCO mapped Jupiter more than 70 times. Each set of 
observations was independently fit. Variations from year 
to year are due to different alignments of the parabola and 
feeds. 

4.2. Optical components 

The primary mirror is a simple offset parabola. In 
the notation of Wollack et al. (1993), the offset angle is 
= 48°, and the focal length is 0.9 m. The dewar views 
the parabola at an angle of = 7° thus the zenith angle 
of the optical axis is 90° - *p -|- *d = 49°. The diameter 
of the parabola in projection is 0.86 m. The rms surface 
tolerance is w 15 yum. 

The chopper design follows that described in Wollack 
(1997) with a few modifications. It is smaller and lighter 
than the Saskatoon chopper, measuring 1.8 by 1.2 m, and 
is driven at the resonant frequency of a flat steel spring 
plate (85 cm by 8.62 cm by 0.167 cm for TOC098) at- 
tached at its middle to the chopper mount. The surface 
tolerance of the plate is w 30 /im. The thickness of the 
spring plate is tuned for each campaign. The resonant 
system produces a sinusoidal sweeping pattern in azimuth 
that requires minimal drive power and produces minimal 
vibrations on the mount. There is no reaction bar on the 
QMAP/TOCO chopper. 

The chopper position is sensed 80 times per chop with 
a set of redundant LVDTs (Linear- Variable-Differential- 
Transformers) , one on each side of the mirror. Errors in 
the sinusoidal chopper motion due to wind loading or in- 
strument ghtches are measured. The temperature of the 
chopper is monitored at all times using nine thermometers 
placed at various locations on the flat. For TOCO, when 
the rms chopper position over a single chop deviates from 
the expected position over the chop by more that 0°.015 
the (lata are rejected. The chopper zero point position in 
azimuth was measured as a function of time in the field 
and found to change less than 0°.01. For QMAP, no data 
cuts were made based on the chopper motion. 

5. OFFSETS AND SIDELOBES 



Emission from the instrument produces signals that can 
potentially complicate the measurement^^ . We generically 
call these signals "offsets" as they are fairly constant over 
long periods and would be present in the absence of a ce- 
lestial signal. Examples are shown in Figure 6. Outside of 
the data selection based on weather, most of the analysis 
effort goes into ensuring that offsets do not contaminate 
the final results. This section addresses the known offsets 
and sets upper limits on their magnitudes for both the case 
of when the source is chopper-modulated (radiation which 
enters the detectors after being affected by the chopper) 
and unmodulated (radiation which enters the detectors di- 
rectly). The offsets from modulated emission can, in prin- 
ciple, occur at any harmonic of the chopper frequency^^. 
In practice, they are predominantly at the lowest spatial 
frequency harmonics. 

In the following we focus on TOCO because the long 
observing campaigns required careful monitoring of the 
offsets. In QMAP, because of the short durations of the 
flights, the offsets were stable. In the mapmaking analy- 
sis, the offsets were projected out of the CMB data using 
a technique described in de Oliveira-Costa et al. (1998). 

5.1. Earth Emission Offset 

Radiation from the Earth can diffract over the front edge 
of the ground screen and enter the receivers. The temper- 
ature of the ground as seen by the detectors is 



Ta 



9ieed{d) 



D = 



An 



1 



± 



with 



cos((7-a)/2) sin((7 -f- Q!)/2) 



(3) 



(4) 



where Te is the physical temperature of the Earth (« 273 
K); .gfood(^) is the gain of the feed as defined by GmaxPn{0) 
where Gmax is the forward gain and Pn{9) is the normal- 
ized beam pattern; fl^ is the solid angle of the Earth 
subtended by the telescope rim; r is the distance from 
the horn to the top front edge of the baffle; 7 and a 
are the diffraction angles; and D is the diffraction coef- 
ficient (Keller 1962). The positive sign in D is for the 
E-field perpendicular to the edge and the negative sign is 
for the E-field parallel to the edge. 

From integrating equation 4 over the geometry of the 
ground screen and the feed pattern, we find the diffracted 
power in D-band into the feed Ta ~ 25 /iK. The front 
baffle, which is in the far held of the feed, has a "Keller 
flare" that reduces the diffracted power over that from the 
sharp edge we assumed in the calculation. Based on our 
experience, calculations of this type involving complicated 
geometries are accurate to roughly a factor of three. We 
have included this factor in Table 6 where results are sum- 
marized. 

We estimate the modulated contribution the same way 
as above but with .gfood(^^) replaced by .gboam(^^), the far- 
fleld gain of the main beam rather than the gain of the 
horn. The calculated power diffracted into the chopper 
sidelobes is « 5 /xK. The front baffle is in the near field 

The Sun is another source. For TOCO, the Sun travels overhead and so camiot be completely blocked at all times. However, for CMB work 
we use data only from when the Sun is fully blocked from the optics. QMAP flew at night. Lunar emission was not seen in either experiment. 
We use the term "harmonic" to refer to the spatial frequencies of the scan pattern as discussed in Section 10. 
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Table 5 
Beam Parameters 



Campaign 


Kal 


Ka2 


Ql 


Q2 


Q3 


Q4 


Dl 


D2 


Predicted 


















O . aCi n— 4 ™\ 


2.76 


2.76 


1.53 


1.53 


1.69 


1.69 


0.124 




^FWHM (deg) 


0.905 


0.905 


0.663 


0.663 


0.702 


0.702 


0.190 




^I'VifM'' (deg) 
Azimuth'^ (deg) 


0.888 


0.888 


0.661 


.661 


0.683 


0.683 


0.192 




203.13 


203.13 


206.75 


206.75 


206.70 


206.70 


205.00 




Elevation (deg) 


41.8 


41.8 


41.8 


41.8 


39.2 


39.2 


40.66 




Polarization 


< — > 


t 


< — > 


t 


< — > 


t 


< — > 




rV\ Cj 1 [CiD ) 


-zi 


-zi 


1 n 
-iy 


1 n 

-iy 


on 


on 
-zU 


oo 
-zz 




Chop ET (dB) 


-47 


-47 


-48 


-48 


-49 


-49 


-54 




QMAP96a 


















Qa (10^'* sr) 


2.83 


2.83 


1.58 


1.58 


1.62 


1.62 








0.931 


0.931 


0.694 


0.694 


0.700 


0.700 






^FiTifM (deg) 


0.882 


0.882 


0.658 


0.658 


0.668 


0.668 






Cross-El (deg) 


-1.35 


-1.35 


+1.35 


+1.35 


+1.35 


+1.35 






QMAP96b 


















Ha (10-* sr) 


2.67 


2.67 


1.43 


1.43 


1.75 


1.75 






^FlVJfM (deg) 


0.932 


0.932 


0.674 


0.674 


0.730 


0.730 






Ofwhm (deg) 


0.831 


0.831 


0.616 


0.616 


0.694 


0.694 






Cross-El (deg) 


-1.35 


-1.35 


-1-1.35 


+1.35 


+1.35 


+1.35 






Elevation (deg) 


41.36 


41.36 


41.37 


41.37 


38.79 


38.79 






TOC097<= 


















Qa (lO-* sr) 


2.76 


2.73 


1.63 


1.74 


1.76 


1.78 


0.183 


0.323 




6% 


6.5% 


4% 


8% 


4% 


6.0% 


5% 


6% 


^fVj/m (deg) 


0.881 


0.871 


0.711 


0.744 


0.718 


0.716 


0.225 


0.306 


^'fVhm (deg) 
Azimuth (deg) 


0.909 


0.909 


0.664 


0.676 


0.711 


0.721 


0.236 


0.306 


203.13 


203.13 


206.75 


206.75 


206.70 


206.70 


205.00 




Elevation (deg) 


41.75 


41.75 


41.85 


41.85 


39.25 


39.25 


40.44 


39.93 


TOC098 


















Ha (10-" sr) 


3.00 


3.00 


1.52 


1.60 


1.76 


1.80 


0.136*' 


0.292 




8% 


8% 


9.8% 


18% 


8.4% 


10% 


5.5% 


5% 


Ofwhm (deg) 


0.860 


0.914 


0.666 


0.669 


0.692 


0.688 


0.201 


0.293 


^FlVffM (deg) 
Azimuth (deg) 


0.907 


0.918 


0.659 


0.681 


0.732 


0.754 


0.194 


0.286 


205.67 


205.67 


209.16 


209.16 


209.06 


209.06 


207.47 


205.73 


Elevation (deg) 


42.05 


42.05 


42.03 


42.03 


39.48 


39.48 


40.63 


40.13 



"Solid angle of beam. 

''The full width at half maximum in the azimuthal and elevation direction. The beam is not symmetric due to smearing in the azimuth direction 
and due to the placement in the focal plane. For QMAP, wc give the major and minor axes of the best fit ellipsoidal Gaussian. 

"^The azimuth and elevation of the beams. For QMAP, we give the relative position in the focal plane. For TOC098, the pointing was slightly 
different. 

■^The edge taper is the ratio of radiant power in the center of the optic to that on the edge expressed in dB. This is determined from a full 

calculation of the current distribution on each optical element using the DADRA program (YRS Assoc.). 
''The chopper was tilted around its normal for TOC097 and so the beam moved in elevation by 0.116 sin(0Q^_j,/jop) deg. There was no tilt in 

TOC098. 

■'^The prediction is for the 1998 configuration. D band is also the most sensitive to slight changes in alignment. 
^The rms in the solid angle during the campaign rounded to the nearest 0.5%. 
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Raided Uiocz 
^ Q3&4 




Fig. 5. — Beam map of focal plane (left) and the synthesized beams (right). Left: Units are degrees on the sky from the center of the focal 
plane. Each beam is normalized with the contours representing 10% in amplitude. The physical separations between the Kal/2 horn and 
the Ql/2 feed is 5 cm. D2 is offset from the focal plane center by 2.9 cm. The distance between the Ql/2 and Q3/4 feeds is also 5 cm. 
This is also a picture of the up/down reflection of the feeds in the focal plane when looking into the dewar. Right: The synthesized beams 
for TOC097 as discussed in Section 10. If the Ka primary beam on the left, for example, is weighted by positive and negative numbers as it 
moves across the sky, one obtains the synthesized Ka pattern on the right. The contours indicate alternate positive and negative lobes. Shown 
are the Ka 9-pt, Q 11-pt, and D 27-pt D beams. From this picture, it is clear that the synthesized beam is sensitive to only a narrow band 
of spatial frequencies. The central dashed line corresponds to RA=0 (Fig. 8) though the synthesized beam location as shown is arbitrary. In 
the full analysis, the beam is smoothed in right ascension. 



Table 6 

Estimated Contributions to the Offsets 



Contributions to the Signal 


D-Band 


Q-Band 


Ka-Band 


Modulated 


Earth emission (directly into feed)^ 


75 nK 


300 ^K 


400 AtK 


No 


Earth emission (in chopper sidelobes)'' 


< 50 fj,K 


< 200 fiK 


< 250 ^K 


No 


Ground screen emission (directly into feed) 


5 K 


6 K 


7 K 


No 


Ground screen emission (in chopper sidelobes) 


80 mK 


80 mK 


80 mK 


No 


Cavity emission, vertically polarized 


3 mK 


13 mK 


15 mK 


No 


Total loading from optics 


5 K 


6 K 


7 K 




Total ground screen emission'^ 


1.5 mK 


1.5 mK 


1.5 mK 


Yes 


Cavity emission (vertically polarized) 


0.3 mK 


1.3 mK 


1.5 mK 


Yes 


Polarized chopper emission'* 


3.5 mK 


1.8 mK 


1.6 mK 


Yes 


(5T|| mirror misalignment in atmosphere 


3 mK 


1.7 mK 


1.7 mK 


Yes 


5T± mirror misalignment in atmosphere 


20 ^K 


10 ^K 


10 fiK 


Yes 


Feed rotation (1°) 


3 mK 


0.8 mK 


0.6 mK 


Yes 


Feed rotation (4°) 


10 mK 


3 mK 


2 mK 


Yes 



"Results of the formal calculation times three to account for modeling errors. 
''Includes a factor to account for near field effects. 

•^Due to a hypothetical temperature gradient across the ground screen. 

■^Magnitude of the largest component of polarized emission. The bold "Yes" indicates that this quantity is computed to 10% accuracy. 
''We estimate that the effective area of the cavity between the chopper and the outer feed baffle is modulated by 10%. 
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of the main beam and so the true values could be up to 
an order of magnitude larger. If the temperature of the 
ground on either side of the telescope were to differ by 10 
K then the offset produced would be < 1 ^K. Similarly, 
any modulation of diffracted power from variations along 
the top of the ground screen are small. From these esti- 
mates we conclude that emission from the Earth does not 
contribute to the microwave signal. 

5.2. Ground Screen Emission Offset 

The antenna temperature from emission of the baffles 

as seen by the receiver is 

Ta « ^^^T^^eba/fieTbaffle^^baffle ~ 5K. (5) 

47r 

where ehaffie is the emissivity of the baffle and Tbaffle is 
its physical temperature. Along with the atmosphere, the 
direct emission from the ground screens is the largest ra- 
diation loading. As with the diffracted groimd emission, 
we ask what portion of this signal is modulated. We find 
that the temperature of the baffle in the main beam, af- 
ter reflecting off the chopper, is « 80 mK. A temperature 
differential of 10 K on cither side of the baffle would then 
translate to an observed offset of « 1.5 mK in the lowest 
harmonics. It is also possible to get modulated emission 
because the polarized emission is a function of the angle 
of the beam with respect to the ground screen. The ge- 
ometry of the enclosure is complicated, but the angle of 
the ground screen implies that the emission will be great- 
est in the vertical polarization. However, the modulation 
will be greatest for the horizontal component. This term 
is similar in character to the polarized emission from the 
chopper but is an order of magnitude smaller. 

5.3. Cavity Emission Offset 

The cavity behind the chopping mirror (Fig. 2)and outer 
feed baffle is effectively black. Radiation from this cav- 
ity reaches the receiver: a) by traveling over the front 
edge of the outer feed baffle and back into the feeds or 
b) through shallow angle diffraction over the outer feed 
baffle onto the parabola and then reflecting back to the 
feeds. This emission is modulated as the chopping mirror 
sweeps back and forth, changing the size and shape of the 
opening to the emitting cavity. Improved shielding of this 
cavity between the two QMAP flights reduced the offset 
at a given angle of the chopper by 1.8 mK or 50% in Ka 
band (Herbig et al. 1998). 

The contribution from radiation traveling over the feed 
baffle and back into the feeds may bo estimated with equa- 
tion 4. For D-band we find that Ta « 3 mK for vertically 
polarized radiation and Ta ~ 0.1 mK for horizontally po- 
larized radiation. 

The contribution from the shallow angle diffraction is 
difficult to compute to even order-of-magnitude accuracy 
because it is so critically dependent on the geometry. We 
note, though, that this term can be of millikelvin magni- 
tude and horizontally polarized. 

5.4. Chopper Polarized Emission Offset 

The emission from the chopping plate is po- 
larized and depends on the angle of the mir- 
ror (Landau & Lifshitz 1960; WoUack et al. 1993 ; 



Cortiglioni 1994). Therefore, as the chopping mirror scans 
across the sky the plate emissivity, as viewed by a feed, 
changes with the position of the chopper. The parallel 
and perpendicular emissivities are: 



eo/cos(^i) & e_L Ri eocos(^i), 



(6) 



where eo is the emissivity at normal incidence and the in- 
cident angle of radiation is 



i(^c) 



COS 



-ki ■ n{(j)c)], 



(7) 



where ki is the propagation vector for the incoming ray, 
n is the normal to the plate, and (j)c is the chopper angle 
in the azimuthal direction. The brightness temperature of 
the plate is 



Temit ~ T'chop [e II I (^'i) I'' 



(8) 



where E^^ and E± are the parallel and perpendicular elec- 
tric field projections onto the normal to the chopping mir- 
ror and Tchop is the physical temperature of the chop- 
per. For the D-band, the peak-to-peak calculated offset 
is Ri 3.5 mK. The dashed curves in Figure 6 show the po- 
larized emission offset calculated for each of the channels. 

In the beam synthesis, the net signal for each chopper 
position is multiplied by the corresponding clement in the 
synthesis vector (Section 10). The offset, since it is an ad- 
ditive signal, is multiplied by the same vectors. It is evi- 
dent from Figure 6 that the offset in the synthesized signal, 
or "synthesized offset," will be larger for the smaller har- 
monics. The polarization directions were chosen to mini- 
mize the synthesized offset for the 3-pt beam as no science 
data is expected from 2-pt beam. 

5.5. Atmospheric Offsets 

When the telescope is properly aligned, the chopper 
scans horizontally through the atmosphere and the atmo- 
spheric emission temperature at all portions of the chop 
is the same. If either the chopping mirror or the entire 
base is misaligned with respect to the horizontal, offsets 
are produced (WoUack et al. 1997) according to 



< Tatm > +- 



'-atm — -'-atm 



5ip\\ -h (5V'_L 



where 



= T,tB.n{e,)sec{e^) 



(9) 
(10) 



is the gradient in the atmospheric temperature, is the 
zenith temperature, and 6z is the fiducial zenith angle of 
the beam. If the chopper is misaligned due to a rotation 
about an axis parallel to the chopper normal, the sky sig- 
nal is changed by 

« 2(50||sin(6i^)sin((?!.c) (U) 

where is the azimiithal chopper angle. If the chopper is 
misaligned due to a rotation about an axis perpendicular 
to the chopper normal, the sky signal is changed by 



-260 I cos(^^2)cos( 



(12) 



Measurements of chopper and base tilt put a limit on the 

L of < 0.1°. Asi 
3 mK and 6Tj_ « 20/iK for 



measured value of S0\\ and 69 ± of < 0.1°. Assuming this 



value, we would expect 6T^ 
the D-band. 
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FlG. 6. — Top Panels: Measured chopper offsets, Oi, in D-baiid as a function of chopper position for a typical night of CMB data. This plot 
shows the raw data (with cal pulses subtracted) synchronously coadded with the chopper azimuth, (pc- Dashed lines arc the calculated offsets 
due to polarized emission from the chopping mirror with the offset at the center of the chop removed, T{4>c)— < T(<f)c = 0) >. The straight 
lines are the calculated offset due to a mirror tilt of S6^^ = 0°.l in the plane of the stationary chopper. When these offsets are multiplied by a 
synthesis vector, the result is constant in time and independent of chopper position. Dl is horizontally polarized; D2 is vertically polarized. 
Bottom Panels: The residual measured offset after the subtraction of the predicted offsets. The form of the residual indicates that the source 
is polarized. The asymmetry in Dl may be due to feed rotation. 



5.6. Observed Offset 

Figure 6 shows the choppcr-position-dependent offsets 
observed in each of the D-band channels on a typical night 
analyzed for CMB observations in the 1998 campaign. 
Plotted in the top frames are the observed offsets (in an- 
tenna temperature) as a function of azimuthal chopper 
angle. Overplotted are the computed polarized emission 
offsets for each of the channels (dashed line) and offset due 
to mirror tilt about an axis parallel to the chopper normal 
(solid line) . The bottom frames show the offsets corrected 
for these two effects. 

There is a clear asymmetry in the offset about (j)c = 
in the top panels of Figure 6. This is most likely caused 
by a misalignment of the chopper with (J^y = 0.1° (equa- 
tion 11), corresponding to the D-band beam centroid mov- 
ing up and down vertically J^'n = 0?0053 as it scans the 
azimuth. This angle is just below the detection threshold 
of our measurements (Section 7). With a zenith temper- 
ature of 10 K in D-band, d9\\ = 0.1° produces a 1.6 mK 
modulation. 

Another mechanism for producing an asymmetry of the 
same magnitude is the rotation of the polarization direc- 
tion of a feed. A « 1° rotation results in a signal of 3 mK 
in D-band (Table 6). We cannot rule out that some part 
of the asymmetry is due to this, though it would be coin- 
cidental to have the asymmetry so similar in both feeds. 
In the HEMT channels, the signature of the asymmetry 
due to feed rotation is opposite in the two polarizations. 
As the chopper sweeps, the emission for one polarization 
goes up while that from the other polarization goes down. 
We see no evidence for such a signal. 

The offsets for Ka and Q band are similar to those shown 
in Fig. 6 and in Fig. 1 of Herbig et al. (1998). Be- 
fore accounting for the chopper emission, the magnitude 
is between 4 to 8 mK. The offsets in Ka and Q bands in 
the Saskatoon experiment were 1/2 to 2 mK, considerably 
smaller than those observed here. We attribute the dif- 
ference to the fact that the TOCO/QMAP system is, by 
necessity, much more compact: the chopper is closer to the 
feed horns and the ground screens are closer to the main 
beam. 



After accounting for the polarized chopper emission and 

the alignment of the chopper, both of which can be com- 
puted accurately, there is still a residual polarized offset 
in all channels of magnitude 2 to 8 mK. In particular, the 
observed residual offset is always largest in the horizon- 
tal polarization as shown for D-band in the bottom panels 
of Fig. 6. The offsets are largest in Q-band, and thus no 
particular characteristic spectrum is identified. This sug- 
gests that the source is modulated cavity emission, thermal 
emission from the enclosure, or a combination of both. We 
have not been able to definitively identify the mechanism 
responsible for the effect though the most likely source is 
the shallow angle diffraction of the chopper cavity emis- 
sion into the parabola. Other potential sources are either 
too small or do not have the correct polarization. 

The offsets we discuss in this section are not atypical for 
CMB experiments. They correspond to the raw detector 
output before any of the symmetries or modulations of the 
experimental design have been utilized. In Sections 11 and 
12 we discuss the offsets after the strong spatio-temporal 
filter of the experimental method has been applied. 

6. ELECTRONICS AND DATA ACQUISITION 

Because it was balloon borne, QMAP was by necessity 
self contained. It had its own command telemetry and 
CCD-based pointing system (Devlin et al. 1998) though 
used the National Scientific Balloon Facility (NSBF) trans- 
mitters to relay data. To operate remotely in Chile as the 
TOCO experiment, a transmitter was added to telemeter a 
compressed version of the data from the telescope on Cerro 
Toco to the ground station in San Pedro de Atacama, 35 
km away. A block diagram of the telescope electronics, 
data system, and telemetry is shown in Figure 7. 

6.1. Electronics 

The eight radiometry signals are processed in the back- 
pack. Each signal is square-law detected with detector 
diodes, amplified, and sent with differential line drivers 
to a processing board. On the processing board, the sig- 
nal is high-pass filtered with a 2-pole RC network with 
/sdB = 1 Hz to remove the DC level which is separately 
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Fig. 7. — Layout of electronics and telemetry system for TOCO. The components on the left are on Cerro Toco and the ones on the right 
are at the ground station in San Pedro (a hotel room in La Casa de Don Tomas), at a comfortable living altitude. The command up/down 
link is a 40 W MotoroUa GM300 radio communicating via a Paccomm packet modem. 



recorded. This filter introduces a small but stable elec- 
tronic offset. The high-passed analog signal is then dig- 
itized using an 18-bit EA analog to digital converter 
The serial digital signal is sent through a shift register and 
latch that converts it to a 32-bit parallel word accessible 
by a computer. The use of SA ADCs is advantageous 
over sampling ADCs because of their superior differential 
nonlinearity specification, which is important when signals 
comparable to or smaller than the 1-bit level are of inter- 
est. Additionally, these ADCs incorporate a digital anti- 
aliasing filter, eliminating any temperature dependence or 
drift in this component of the electronic band pass. 

In the CMB analysis, we use HEMT data from 14 to 
55 Hz ( 4 < ; < 200 ) and SIS data from 15 to 60 Hz ( 
4 < I < 500 ). The electronic bandpasses are defined on 
the lower end by the high-pass filter and on the high end 
by the EA chip. Over our frequency range, the phase re- 
sponse of the chip is linear with frequency and so the EA 
introduces a time delay in the signal. The amplitude re- 
sponse over the CMB frequencies varies by less than 0.4% 
over this range. 

6.2. Computers and Telemetry 

Two single board computers, which handle the data and 
the pointing, are located on the telescope. They are syn- 
chronized with a common clock and communicate remotely 
with the computers operating at the base station. In the 
following, we focus on the configuration for TOCO. 

The "Data PC" logs the detector output and the po- 
sition of the chopping mirror as well as various voltages, 
currents, temperatures, etc. which we use to monitor the 
telescope. Two versions of the data are recorded. A com- 
plete version (w 1 G byte/day), which is used in the final 
analysis, is stored on a 4 G bytes hard drive in the data 
computer. This drive is contained in a pressurized vessel 
to prevent damage resulting from operation at high alti- 
tude These drives are swapped out every two to four 
days. The data are then uploaded onto the computer sys- 
tem in the ground station, and stored on Exabyte tapes. 
A second compressed version of the radiometry and house- 
keeping data is bi-phase encoded and sent real time to the 
ground station. 

The "Pointing PC" records the position of the 17-bit ab- 

23 One half of an Analog Devices AD1878 

■^'^ When not enclosed in a pressurized and dust-free container, most hard drives, especially high capacity ones, were found to fail on a timescale 
of a day. 

■^^ We decided against observing on both sides of the South Celestial Pole (which would have produced interlocking scans thereby facilitated 
map production) to maximize the stability of the instrument and to minimize the complexities of the analysis. 



solute digital encoder on the telescope azimuthal bearing, 
controls and monitors the azimuthal drive motor, records 
the time from a GPS receiver, and interprets commands 
sent remotely from the ground station. A CCD camera 
and Matrox digital image processing board can track stars 
during calibration and pointing verification procedures. 
The command status and all other information on the 
pointing PC is passed to the data PC for logging. The 
Pointing PC can also reboot the Data PC, an operation 
we sometimes found necessary. 

Two radio links allow us to communicate with the tele- 
scope from the ground station. A high frequency link at 
1.4 GHz (bandwidth of 100 kHz) with a 2 W transmitter 
links the data computer with a ground station PC provid- 
ing the bi-phase data. A marine radio operating at 460.5 
MHz communicates via a packet modem with the pointing 
computer, providing commanding. 

The ground station computers receive and store the 
data, archive the data to tape, and run the commanding, 
display, and alarm software. From the ground station, the 
telescope can be slewed in azimuth, chopper parameters 
can be adjusted, and the cooling power (to stabilize the 
temperature of the warm electronics) can be increased or 
decreased. Most major systems can be turned on and off 
remotely. 



7. OBSERVATIONS 

While both QMAP and TOCO were designed to mea- 
sure the anisotropy, their approaches were completely dif- 
ferent. QMAP was designed to make a true map of the 
sky. The data from the time stream were pixelized on the 
sky and the analysis was done on the sky map. In the first 
flight of QMAP, the chopper swept horizontally at 4.7 Hz 
and the gondola wobbled in azimuth with a period of 100 s 
about a meridian containing the North Celestial Pole. This 
gondola motion, combined with the rotation of the Earth, 
produced a highly interlocking scan pattern that allowed 
for the clean separation of instrumental effects from the 
celestial signal (de Oliveira-Costa et al. 1998). 

The TOCO scan was designed to measure the angular 
power spectrum^^. The telescope optical axis is fixed in 
azimuth and elevation, as indicated in Tables 1 & 7, and 
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Table 7 

Telescope and Chopper Motion 



Parameter 



QMAP96a 



QMAP96b 



TOC097 



TOC098 



Chopper frequency, fc (Hz) 4.7 4.6 4.6 

Chopper amphtudc in azimuth, fflc(max) (deg) ±10 ±2.5 ±2.96 

Elevation at center of chop, del (deg) 40.7 40.1 40.5 

Azimuth of center of chop f 6az (deg) scanned scanned 204.9 

Amplitude of wobble (s) 100 50 none 

Amplitude of wobble in azimuth (deg) ±5° ±1?5 none 



3.7 

±2.02 

40.63 

207.47 

none 

none 



"This is the physical motion in a horizontal plane. The amplitude on the sky is 2(pc cos 9f,l ■ 



the chopper sweeps the beam across the sky. The beams 
cover an annulus around the SCP as shown in Figure 8. 

The TOCO observing site is located at an altitude of 
5200 m on Cerro Toco in the Northern Atacama desert 
in Chile near the borders of Argentina and Bolivia. The 
latitude is 22°.95 South and the longitude is 67°. 775 West. 
A building for an abandoned sulphcr mine blocks the oc- 
casional 60 kt winds. The Atacama is one of the highest, 
driest deserts in the world and is therefore a good place 
for millimeter and centimeter wave observations. We find 
that the weather is good enough for D-band CMB obser- 
vations Ri 50% of calendar time between September and 
January. 

8. CALIBRATION 

Calibration uncertainty is the largest source of exper- 
imental error for both QMAP and TOCO as it is for 
many other CMB experiments. In general, calibration is 
dominated by systematic effects. Multiple independent 
measurements are required to determine the uncertainty. 
Short of this, one's knowledge of the experiment is called 
upon to set the uncertainty and its distribution 

A calibration for TOCO involves the following steps: a) 
position the telescope so that Jupiter either rises or sets 
through the azimuthally swept beam; b) make a map of the 
source and find the best fit amplitude, position, and beam 
profile; c) compute the brightness of Jupiter accounting for 
its position with respect to Earth; d) correct for slow drifts 
in calibration, correct for the electronic response of the 
receivers, and convert Rayleigh-Jeans to thermodynamic 
temperature. A similar procedure was used for QMAP. In 
the following, we consider each element of the calibration 
and its associated uncertainty. 

There are two general classes of calibrators, those 
that fill the beam and point sources. FIRS 
(Ganga et al. 1993 ), BOOMERanG, and MAXIMA 
used the dipole as calibrated by both FIRAS and DMR 
aboard COBE. The dipole signal corresponds to angular 
scales larger than, and detection frequencies lower than, 
those used for CMB data analysis. Thus knowledge of 

^® One may, for example, conservatively use the diflference between two measurements as the uncertainty (e.g., de Bernardis et al. 2000) even 
though the formal uncertainty is the difference divided by \/2. To determine error bars on quantities derived from the data, one marginal- 
izes over the distribution that describes the calibration. Generally the distribution is taken as Gaussian, though one might reach different 
conclusions if a different distribution in fact described the data. 
27 3C461, 2321-|-583IAU(1950), I = 111?7, b= -2?1 

2^ A typical frequency dependent loss in the system has a negligible ailect on the calibration. 



the electronic transfer function, beam profile, and spa- 
tial filtering are necessary for the extrapolation. At the 
other extreme, point sources are at smaller angular scales 
and higher post-detection frequencies than those used for 
CMB analyses. They have the advantage that the calibra- 
tion and beam profile are measured simultaneously. Still, 
though, one must account for the electronic transfer func- 
tion and spatial filtering. 

QMAP was calibrated with Cas-A^^ and TOCO was 
calibrated using Jupiter. As the beam sweeps across the 
source, the detector output is given by 

CV{t) = l-[ [ [ A,{v)g{v)S,{VL)Pn{v,x-^{t'))x 

^ Jv Jfl Jt' 

I{t - t')dxdvdt' (13) 

where V is the voltage out of the detector, Ae is the ef- 
fective area of the telescope dish, g{u) is the passband of 
the receiver, Sv is the source surface brightness in units of 
W/m^ srHz, P„ is the normalized gain of the antenna, x 
is the direction on the sky, and / is the impulse response 
function of the electronics^*. Here C is the calibration con- 
stant that relates the source temperature to the measured 
voltage. We do not explicitly account for sampling in the 
integral. 

Each term in equation 13 has an error associated with 
it that contributes to the net uncertainty. The effects in- 
clude: a) (Jsi the intrinsic uncertainty in the brightness 
temperature of the source at the calibration frequency ex- 
trapolated from previous measurements; b) ay, the uncer- 
tainty in the measured temperature of the source; c) ctq, 
the Tinccrtainty in the measured solid angle of the beam 
wliicli includes any beam smearing; d)the uncertainty in 
the measured receiver bandpass, Ug, and center frequency, 
<yv', e) (J drift, the uncertainty in the change of the calibra- 
tion between when it is measured and when it is applied; 
and f) (7/, the uncertainty in the frequency response of 
the instrument. Items (b) and (c) are derived from mea- 
surements of the source and are covariant; generally one 
conservatively treats them as independent errors. Com- 
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Table 8 
Observing Parameters 



Campaign 


Kal 


Ka2 


Ql 


Q2 


Q3 


Q4 


Dl 


D2 


OMAPQfia 


















Sample frequency, fs (Hz) 


752 


752 


752 


752 


752 


752 






Samples per chop, Ns 


160 


160 


160 


160 


160 


160 






Ofwhm per scan 


20 


20 


30 


30 


30 


30 






Samples per Ofwhm (F/S)^'*' 


0.78/3.2 


.78/3.2 


0.56/2.7 


0.56/2.7 


0.58/2.7 


0.58/2.7 






Point Source Sensitivity T{uc){fJ.KJy~^) 


128 


128 


132 


132 


129 


129 






v^iViAi yoD 


















SaTTiTilp frpnnpncv f 


736 


736 


736 


736 


736 


736 






Samples per chop, Ns 


160 


160 


160 


160 


160 


160 






Ofwhm per scan 


5 


5 


8 


8 


8 


8 






Samples per Ofwhm (F/S) 


3.1/6.5 


3.1/6.5 


2.4/5.7 


2.4/5.7 


2.3/5.6 


2.3/5.6 






Point Source Sensitivity T(i'c)(tJKJY~^) 


136 


136 


146 


146 


120 


120 






1 uuuy I 


















Sample frequency, j s (Hz) 


oDo 


oDo 


ODo 


oDo 


ODo 


ODo 


1 A 70 




Sky detection frequency (Hz) 


lo-OO 


lo-OO 


lo-D4 




Z0-D4 


Z0-D4 






Samples per chop, Ns 


CO 

oU 


en 


en 
oU 


oU 


QO 
oU 




QOO 

ozU 




Ofwhm per scan 


6 


6 


9 


9 


9 


9 


30 




Samples per Ofwhm''(F/S) 


1.3/3.0 


1.3/3.0 


1.0/2.6 


1.0/2.6 


1.0/2.6 


1.0/2.6 


1.1/5.2 


1.7/6.7 


Point Source Sensitivity r(!/c)(A'KJy~^) 


131 


131 


128 


120 


119 


118 


87.3 


49.5 


Tocogs" 


















Sample frequency, js (Hz) 


296 


296 


296 


296 


296 


296 


1184 


1184 


Sky detection frequency (Hz) 














19-59 


19-63 


Samples per chop, Ns 


80 


80 


80 


80 


80 


80 


320 


320 


Ofwhm per scan 


4 


4 


6 


6 


6 


6 


20 


20 


Samples per Ofwhm (F/S) 


1.9/3.6 


1.9/3.6 


1.5/3.2 


1.5/3.2 


1.4/3.0 


1.4/3.0 


1.7/6.7 


2.5/8.1 


Primary Aperture Efficiency, r\p 


0.58 


0.58 


0.65 


0.62 


0.56 


0.55 


0.52 


0.24 


Chopper Aperture Efficiency, r/p 


0.10 


0.10 


0.11 


0.11 


0.10 


0.10 


0.087 


0.040 


Point Source Sensitivity r(i/c)(AtKJy~^) 


121 


121 


138 


131 


119 


116 


117 


54.7 



"^F/S refer to the fast and slow parts of the sinusoidal motion near the center and edges of the chop. 

'The sky is not Nyquist sampled at the center of the chop. This is accounted for in the mapmaking and beam synthesis. The point of the 
beams, = \/161n2/(9^vi'£rA/i is = 212,273, and 955 for Ka, Q, and D respectively. The undersampling for QMAP96a is severe and limits 
the map reconstruction in the current pipeline to I 180. 




Gahtxy Cu! 

Fig. 8. — Scan pattern for the TOC098 campaign around the SCP. The scans are shown as straight lines for each feed for 15 minute 
intervals around the sky. A more detailed picture would show the lines with slight curvature to account for the projection. The center of the 
chop is at an azimuth of 207°. 5 and an elevation of 40°. 6. The chopper sweeps out 6? 12 on the sky as it scans in azimuth. This results in a 
ring around the SCP approximately centered on 5 = —61°. The map of the galaxy is from Schlegel et al. (1997). 
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bining these sources, the intrinsic cahbration is given by: 




We evaluate these terms in the following and summarize 

the results in Table 9. 

8.1. Calibration Source Brightness Temperatures 

The flux from a source is typically given as f^, = 
J Sv{9,(p)dQ, and is reported in units of Janskies 
(1 Jy = 10-26 W/m^Hz) function of frequency. A 

power-law model is used to extrapolate the calibration fre- 
quencies to the observation frequency, = fo{i//i^o))^. 
When the source angular size is a significant fraction of 
the beam or is changing in time, as it is for planets, it is 
more convenient to use the brightness temperature. 

8.1.1. Cas-A Flux for QMAP 

Cas-A is unresolved at the QMAP beam size 
and the correction for its finite size is negligible. 
From a combination of the data (Baars et al. 1977; 
Chini et al. 1984; Mezgcr 1986), we find log(5^/Jy) = 
(5.713±0.023) - (0.759±0.006) log(i//MHz) at epoch 1980 
(Herbig et al. 1998). At 31.3 GHz, the mean frequency of 
the QMAP Ka bands, the flux corresponds to 199.9 Jy. 
Baars et al. (1977) give the percentage annual decrease 
as SS^/S^ = 0.97(±0.04) - 0.3(±0.04) log(i^/GHz) and so 
this value is reduced to 183 Jy for epoch 1997. When one 
takes into account all of the above errors, including cr^, the 
uncertainty is {(Js/S) = 8.7% for both Ka and Q bands. 

After the QMAP data release, Mason et al. (2000) re- 
ported Scas,w%s, — 194 ± 4.7 Jy at 32 GHz which we con- 
vert to 195 Jy for epoch 1997. This measurement, which 
is very close to our frequency, greatly reduces the uncer- 
tainty associated with the interpolation and secular de- 
crease. When the errors in the central frequency are in- 
cluded (Table 1), the uncertainty is [as IS) = 2.7%. For 
Q band, the flux is 159 ± 4.8 Jy. The slightly larger error 
is from the extrapolation from 32 to 41 GHz. The result 
is an increase in the temperature scale of the QMAP data 
by 6.6% and reduction in the calibration uncertainty, as- 

8.1.2. Jupiter Temperature for TOCO 

The brightness temperature of Jupiter is measured by 
Ulich (1981), Griflin et al. (1986), Ulich and Mason et 
al. (2000), and is taken to be 152, 160, and 170 K in 
Ka, Q, and D band respectively with an intrinsic calibra- 
tion error oi (Js / S — b%. As the temperature is a weak 
function of frequency across our bands, the uncertainties 
resulting from and ap are negligible. The Jupiter cali- 
bration temperature is obtained by scaling the brightness 
temperature by the ratio of the solid angle of Jupiter (de- 
termined from ephemerides) to the measured solid angle 
of the beam. A typical measured temperature of Jupiter 
is 15 mK, 30 mK, and 350 mK in Ka through Dl bands 
respectively. 

Finally, one applies a correction to convert small changes 
in antenna temperature to small changes in thermody- 
namic CMB temperature. We use: 

( dT^\ x^e- 



where x = hv/kTcMB- For the TOC098 D-band data, 
the data are multiplied by r]~^, or 1.66 for Dl and D2 to 
convert from data calibrated on Jupiter's brightness tem- 
perature to thermodynamic units referenced to the CMB. 
For the Ka and Q band data we multiply by 1.02 and 1.05 
respectively. The error in these measurements depends on 
the knowledge of the band centroids and introduces < 1% 
uncertainty. 



8.2. Measured Beam Solid Angle and Temperature 

To convert the measured fluxes to a temperature, the 
beam solid angle must be known. Table 5 gives the 
beam determination for all campaigns along with the re- 
sults of a computer model. For QMAP, the solid an- 
gle was determined from one in-flight mapping of Cas-A 
(Herbig et al. 1998). The statistical error on the fit var- 
ied between 1% and 3% depending on the fiight and band. 
As the beam fitting includes modeling of the instrumen- 
tal offsets, there is additional systematic error resulting 
in a net uncertainty of trn/n = 5 %. For TOCO, the 
beam solid angles were determined from a global fit to the 
Jupiter calibrations made during good weather. This was 
done separately for both campaigns. In total, Jupiter was 
mapped over 70 times. All results from the four campaigns 
are consistent with our models. From all of our measure- 
ments of the beams, we conclude that the error on the solid 
angle for Ka and Q bands is 5%, for Dl it is 5%, and for 
D2 it is 5.5%. These values are dominated by systematic 
errors. 

It is often convenient to parameterize the beams with a 
two dimensional Gaussian profile. For Ka and Q bands, 
this introduces a negligible error. We tested for this in 
TOC097 (Torbet et al. 1999) analysis where we used the 
measured Ka and Q-beam profiles in place of the param- 
eterized profiles and found < 1% difference in the final 
results. The D-band channels are less well approximated 
by Gaussian profiles and thus one must use the measured 
profile for accurate results as was done in Miller et al. 
(1999). 

For each calibration, a seven parameter model is fit to 
the data. We find the position of the source, the best 2-D 
Gaussian parameterization including orientation, the am- 
plitude, and an offset. For QMAP, the 1-2% statistical 
error on the amplitude is dominated by a 3% systematic 
uncertainty in the algorithm to extract the amplitude. In 
the TOCO experiment, the standard deviation of the fit- 
ted amplitude is 4-10% for all the HEMT channels in both 
seasons. The variance is a result of atmospheric fluctua- 
tions and finite HEMT sensitivity and thus averages down 
as the square root of « 20 independent high quality maps 
in each season to a value of 2%. There is a small additional 
uncertainty due to the fact that not all fitting algorithms 
give the same results. The net result is to increase the 
uncertainty in the amplitude to av /V = 3%. 



8.3. Calibration drift 

The physical temperature of the TOCO instrument can 
vary by 50 K in a day. Even though all critical components 
are thermally regulated, there are temperature changes 
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that lead to changes in gain. Jupiter was observed, on 
average, within two hours of the beginning of the CMB 
observing time. Changes in system gain on time scales 
shorter than 24 hours were monitored with an internal 
calibration signal ("cal pulse") with an effective temper- 
ature of ~ 1 K in all bands. This pulse was turned on 
40 msec twice every 200 seconds. The amplitude of the 
pulse was regressed with the body temperatures of the 
noise sources, warm electronics temperatures, and cryo- 
genic temperatures in all bands. The fit coefficients are 
consistent with laboratory measurements and show that 
the pulse amplitude is constant but that the system gain 
is a function of the microwave amplifier temperatures. 

For TOCO, a typical long term (50 days) variation is 
15% in Ka band, 5% in Q band, and 20% in D-band. 
The cal pulse amplitudes follow the general trends in the 
Jupiter calibrations in all bands. Prom the cal pulses 
and the Jupiter observations, we derive a calibration drift 
model which we apply to the data. In D-band, atmo- 
spheric fluctuations made use of the cal pulses to correct 
drifts over periods of less than six hours problematic. The 
uncertainty in the model is estimated to be 5%. 

A similar approach was taken with QMAP though the 
cal pulses were clearer and so the uncertainty in correcting 
for the 5% drift is negligible. 

8.4. Electronic passband 

Observing a point source is similar to exciting the elec- 
tronics with a pulse. Consequently there are frequency 
components up to / « 'i[(j>cCos{6ei)0FWHM].fc (Table 7) 
of the post-detection electronics. The CMB signal is at 
comparatively lower frequencies. We model the full elec- 
tronic response of the system, including sampling, and find 
that the CMB data should be reduced by 1.7% in D-band, 
1.5% in Q-band, and 1% in Ka-band for TOCO and 1.5% 
for QMAP over what one would get without the correc- 
tion. We estimate the imccrtainty in this to be 1% for 
TOCO and 5% for QMAP. These shifts were not reported 
in the original papers as they were much smaller than the 
uncertainty, though we include them here. 

Prom Table 7, one sees that we are slightly undersam- 
pled during the fastest part of the chop (where the least 
amount of time is spent). To check for a possible sys- 
tematic effect in TOCO associated with this, calibrations 
were done with the source at the center of the chop and 
off to one side and with different chopper amplitudes. The 
results of these tests arc statistically consistent with the 
nominal calibrations. For TOCO, the CMB anisotropy re- 
sults are insensitive to the slight undersampling because 
the CMB detection frequency is far below the sample 
frequency-due to the beam synthesis-and because in the 
beam synthesis wc simulate the sky scan. 

For the first flight of QMAP the undersampling is more 
drastic and was not included in detail in the mapmaking. 
Thus, the QMAP data should only be considered valid up 
to I = 180^'^ For both flights, the calibration data were 
processed in a manner similar to the mapmaking and so 
beam smearing effects were accounted for in an average 
sense. 

The phase response of the full instrument as a function 
■^^ Neither the undersampUng nor the cahbration shift were a<;counted 
Atmospheric fluctuations preclude direct mapping from the ground. 



of frequency was measured in the lab and determined from 
observations of the Galaxy and Jupiter. We find that the 
phase is linear over the range of frequencies applicable to 
both CMB and point source observations. 

8.5. Combining the calibrations between bands and 
systematic shifts from previous results. 

Both QMAP and TOCO have multiple detectors, the 
measurements of which are combined into one angular 
spectrum. The net calibration uncertainty is a combi- 
nation of terms that are completely correlated between 
channels, such as as, crn, cr drift, & c/, and terms that 
are uncorrelated, such as try and cr,,. When the data are 
combined the last two terms become negligible. The uncer- 
tainties are 8% for D1+D2, for the TOCO HEMTs, and for 
the combination of the Dl+D2+HEMTs. For QMAP the 
net uncertainty is 7.6%. In the regions where the QMAP 
and TOCO angular spectra overlap, only as and an are 
correlated and the combined uncertainty is 6.4%. 

These uncertainties are slightly different than those pre- 
viously quoted for these experiments and are the result 
of a complete reassessment of the calibration errors. To 
correct the previously published results, the QMAP data 
should be multiplied by 1.051, The TOCO D-band data 
(Miller et al. 1999) should be multiplied by 0.983, and the 
TOCO HEMT data (Torbet et al. 1999) by 0.99. 

9. OBSERVING THE ANISOTROPY 

The anisotropy is a two dimensional random field in 
temperature. The goal of CMB anisotropy experiments 
is to measure the characteristics of that field. The three 
methods in use are direct mapping, time-domain beam 
synthesis, and interferometry. For a small number of de- 
tectors, direct mapping makes the most efficient use of 
them, beam synthesis is the next most efficient, and inter- 
ferometry is the least efficient, as wc discuss in Section 10. 
The best strategy to use, though, depends as much on the 
control of potential systematic errors as on raw sensitivity. 

The QMAP experiment was designed to make a direct 
map. By this we mean that the time ordered data are as- 
signed a sky pixel number as they come out of the detector. 
Slow drifts in the detector output may be removed from 
the map using a variety of methods (Cottingham 1987; 
de Oliveira-Costa et al. 1998; Hivon et al 2001). For ah 
methods, though, a heavily interlocking scan strategy is re- 
quired for robust results. QMAP produced multiple » 0?8 
resolution maps of the CMB from two flights. These maps 
were found to be consistent with each other and with the 
Saskatoon maps (Xu et al. 2000). In other words, two en- 
tirely separate experiments measured the same temper- 
ature variations in the sky in overlapping regions. The 
other mapping experiments to do this are COBE/DMR 
(Smoot et al. 1992) with FIRS (Ganga et al. 1993 ), and 
COBE/FIRAS with COBE/DMR (Fixsen et al. 1997). A 
wide range of systematic checks have been applied to 
QMAP and it passes them all. 

The TOCO experiment used time-domain beam synthe- 
sis^°. To our knowledge, the method was first employed in 
Netterfield et al. (1995). Since the Saskatoon experiment, 

for in Xu et al. (2001) and Wang et al. (2001). 
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Table 9 

Summary of Contributions to Calibration Uncertainty for Individual Channels 



Channel as (%) an (%) ay (%) adrift (%) a„ (%) ai (%) ac (% ) 

Dl 5 5.5 3 4 1 1 9.1 

D2 5 5 3 4 1 1 8.8 

TOCO HEMTs 5 5 3 4 • • • 1 8.7 

QMAP HEMTs 2.7 5 3 • • • • • • 5 8.1 



The intrinsic Jupiter calibration uncertainty is given by as (Section 8.1), (tq is the uncertainty in measured solid angle (Section 8.2), ay is 
the uncertainty in measured Jupiter brightness temperature (Section 8.2), and adrift (section 8.3) is the uncertainty due to the calibration 
drift between observations, arj refers to the Rayleigh-Jeans to thermodynamic conversion due to the uncertainty in the centroid, and crj is the 
uncertainty on the electronic bandpass correction. These calibration errors are in temperature. 



we have refined the technique, applied it to a multifeed 
system, and incorporated numerous cross checks and sys- 
tematic checks of the robustness of the solution. Most of 
the remainder of this paper is devoted to describing those 
checks. 

10. TIME DOMAIN BEAM SYNTHESIS 

As with all ground-based CMB experiments, the effects 
of atmospheric fluctuations must be strongly suppressed. 
TOCO uses the chopping flat to scan the beam across the 
sky in a sinusoidal pattern (Section 4.2). In a postde- 
tection analysis, the time ordered data are multiplied by 
synthesis vectors, SVn,i, that have half the period of the 
chopper cycle. Thus, the data with the chopper moving in 
one direction is coadded with the data with the chopper 
moving in the other direction. We form 

Nc 

tn = J2^^n,idi (16) 
i=l 

where rf, is the vector containing the raw data from a full 
chop cycle, n is the index or "harmonic" of the synthesized 
beam, and Nc is the number of samples in a chop cycle. 

There is no set prescription for SVn,i. The best choice 
depends on the scan pattern (e.g., sinusoidal or triangu- 
lar), the desired degree of orthogonality between synthe- 
sized beams, the shape of the resulting window function, 
and the orthogonality to any potential offset. For instance, 
one may pick SVn,i so that f„ is insensitive to the secant 
dependence of the atmospheric gradient. For the sinu- 
soidal scan patterns, we flnd a useful set is given by 

SVniOC cos(7r(n-l)(l-sin(27r(i-l/2)/iV))/2), i = l,...Nc 

(17) 

though ultimately we tune the SVn.i- The synthesis vec- 
tors are effectively apodized sine functions each of which is 
designed to produce a different effective antenna pattern, 
and thus probe a different angular scale. 

The phase of the electronic signal with respect to the 
position of the beam on the sky is determined by forming 
the quadrature signal (data from the first half of chopper 
sweep minus data from the second half resulting in mini- 
mal sensitivity to celestial signals) as a function of phase 
for each harmonic. Then the phase that nulls the galac- 
tic signal over twenty-five of the best observing days is 
found. There is a small harmonic dependent component 
to the best fit phase that is well modeled with a linear 



fit. The phase shift is incorporated into equation 17 for all 
analyses. The galaxy- null determined phases agree with 

the phases determined from Jupiter observations, derived 
with completely independent code. In addition, the entire 
data analysis is redone after setting the phase ahead and 
behind the best fit by twice the error derived in the fit. 
No changes in the final results are seen. 

The resulting effective antenna sensitivity patterns, or 
synthesized beams i?(x), are given by 

H{ic) = (^SViG{ic-Xi)) (18) 

\ j / RA bin 

where the center of the main beam is located on the sky 
at position Xi, SVi is the synthesis vector, and 

is the main beam pattern of the telescope pointed at the 
center of the chop. Here x measures the position in the 

azimuthal direction, y changes with elevation, and the 
beamwidth of the telescope is a = 0FWHM/v'81n2. Beam 
smearing due to the finite size of the sky bins is incor- 
porated into equation 18. As tests, the A/D sampling 
and the change in the effective horizontal beam width as 
a function of chopper speed arc included in equation 19. 
These effects are negligible. As noted earlier, equation 19 
was replaced with the true beam profile for D band. This 
was not necessary for the other bands. The synthesized 
beams are then normalized, by adjusting the amplitude of 
SVn,i, so that: 

J \H{x)\dx = 2 (20) 

In summary, the synthesis vectors and beam synthesis in- 
corporate all known aspects of the motion of the beam on 
the sky. 

Data are binned according to right ascension at the 
center of the chopper sweep into 768 "fundamental bins" 
around the complete circle shown in Figure 8. The funda- 
mental bins are subgrouped into NjiA right ascension bins 
where Njia depends on the harmonic. Thus, the average 
in equation 18 depends on harmonic. To avoid statisti- 
cal bias, it is important that the mean, sample variance, 
and error on the mean in each of the Nra bins is com- 
puted for all the data that land in a bin in a given night. 
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For example, for the 1998 Q-band beam, there are 

■^RA = 20 bins around the circle and so in one night, 
?a 1.7 X 10^ values of t4 (equation 16) are averaged to- 
gether for each bin. At low I, the variance in any RA bin 
is larger than that expected from detector noise alone due 
to atmospheric fluctuations. For integrations longer than 
roughly 3 minutes, the noise is stationary for the cuts de- 
scribed below. Finally, the results from individual nights 
are averaged. 

The harmonics are analyzed individually at low I or 
in groups at higher I. Data from different detectors are 
also combined. The full theory covariance matrix, Ct, is 
computed for each harmonic or combination (Bond 1996) 
along with the Knox filter (Knox 1999) and effective spher- 
ical harmonic index, Ig, of the observations. The noise co- 
variance matrix, Cjv, is computed from the data. Other 
than detector noise, which is uncorrelated between bands, 
the dominant contribution is the atmosphere. Unlike 
Saskatoon, the frequency bands are not subdivided and 
there is no East- West chopping. 

For each group of harmonics, we find the likelihood as 
a function of the band power, ST^, according to^^: 



L{ST^) = 



1. 



27r^/2|M|i/2 



exp(-i^M-H/2), 



(21) 



where M = Ct{ST^) + Cm- All the analysis is done as a 
function of 5Tf, though we report 5Ti because it gives a 
direct measure of signal-to-noise as the detector output is 
proportional to temperature. 

It is sometimes convenient to estimate the signal-to- 
noise for a single measurement for a single harmonic. 
This can then be generalized to a set of N independent 
measurements (Knox 1995). The measured rms ampli- 
tude of the sky fluctuations is given by A = STi^J I(W) 
where I{W) = Y^i{Wi/l), with Wi the window func- 
tion, encodes the efficacy of the synthesized beam. The 

noise of any measurement is given by KnT/r^l^ where 

1/2 

, Tabs is the time spent observing 



a point, and T is given by equation 2. From fits to the com- 
puted I{W), we find ^/I(Wj « 2.6/n°-^5 „ e/a^-^s/O-^s 
where n denotes the n-pt function, at is the Gaussian 
width of the beam, and le is the associated effective spher- 
ical harmonic index. The signal-to-noise for a particular 
measurement is then 



S/N 



We may interpret this result as saying that for a fixed 
sensitivity, T, and a flat angular spectrum, (5T(, beam 
synthesis reduces the effective temperature of the sky by 
^TDBS ^ 0.7/afH°-^. This form fits Ka-band through 
D-band data to 30% accuracy for TOCO 

It is worth contrasting beam synthesis in interferom- 
etry with time domain beam synthesis. In interfer- 
ometry, each baseline yields one synthesized beam en- 

Following the notation in Ncttcrficld et al. (1997), wc call a synthesized beam with n lobes an n-pt beam (see Figure 5). 

In this paper, we report the band power following Bond as opposed to Netterfield. The diflference is 2{l + l)/(2i -|- 1) ~ 0.2 % at i = 200 

For the classic single diflference, k = \/2{l'^ + (—1)^] = 2. 
^'^ Other functional forms work as well. This is presented to aid in estimating the S/N. For SK, which used diflferent criteria to synthesize the 
beams, 'yTOSS _ 3. 6/(70. 5^^^ roughly a factor of two higher at low (. However, the window functions were less well localized. 
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veloped by the primary beam pattern of a single ele- 
ment. A similar expression to equation 22 obtains for 
interferometers with Gaussian main beams of width u^"* 
(Hobson et al. 1995; White et al. 1998). We find 



S/N 



to 



■1 



1/2 



(23) 



where Iq = 27ruo is determined by the separation of the 
two antennae. For the 5-pt Q-band synthesized beam, 
Ofwhm = 0°7 ( ab = 0.0052 ) and le = 87. Similar cov- 
erage in l-spacc would be obtained with an interferometer 
with 0FWHM ^ 4?86 ( fj^"* = 0.036) and = 87. From 
these, 7^-035 ^ o.27 and 7^"* = 0.45. To the level of 
accuracy of the fitting functions, these are equivalent. 

The sensitivity advantage of time domain beam syn- 
thesis is that multiple n-pt functions arc measured si- 
multaneously with a single detector. In the parlance of 
Fourier transform spectroscopy, there is a multiplex or 
"Felgate" advantage over an interferometer with just a 
few antennae. However, as the number of interferometer 
baselines scales as Uaijia — l)/2, where Ha is the number 
of antennae, large interferometers achieve high sensitiv- 
ity (Padin et al. 2000; Pryke et al. 2001). For the ideal 
mapping experiment, with minimal baseline subtraction, 
the advantage over both interferometry and time domain 
beam synthesis is that in a fixed amount of time more 
spatial modes can be measured. 

11. DATA SELECTION 

Most of the analysis effort goes into data selection and 
testing to make sure that the selection does not bias the 
final result. The largest cuts remove data contaminated 
by the galaxy, by the atmosphere, and by unstable off- 
sets. Partial descriptions of the cuts are given in Miller 
et al (1999) and Torbet et al. (1999). In this section we 
describe the cuts for the TOC098 D-band data and the 
consistency checks as they arc indicative of the process for 
all channels. A summary of the cuts is given in Table 10. 

11.1. Cuts to the time line 

The initial cuts are made to the raw time-ordered data 
in order to excise extreme events, such as a nearby ob- 
ject entering the beam or very bad weather. Data are 
examined in 6.5 second averages (24 chop averages with 
fc = 3.7 Hz). The internal calibration pulses are removed 
and a rough cut is made at a nominal rm,s level based on 
the long term observing characteristics. This cut removes 
w 50% of the data. 

The next set of rms cuts is made to the synthesized 
data. In general, the higher harmonics are less sensitive 
to atmosphere and therefore require less severe cuts. The 
data are binned into fifteen minute averages and, for each 
harmonic on each day, the minimum value of the average 
rms is selected as the baseline. For Dl, all data within 
25% of the minimum are accepted for harmonics 8-21. For 
harmonics 5-7 everything within 20% is accepted. For the 
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case of D2, harmonics 8-17 are cut at 30%, and 5-7 are cut 
at 20%. Harmonics < 4 are rejected entirely as they are 
corrupted by atmospheric fluctuations. 

We analyze 28 nights for Dl and 23 for D2. (The num- 
ber is higher for Dl due to high offsets in D2 on several 
nights leading to the decision to exclude these data from 
the analysis.) In order to prevent signal contamination 
from times of large atmospheric fluctuations, the previ- 
ous and subsequent fifteen minute segments are eliminated 
from each segment cut by the above criteria. The effect 
of this cut is to keep 5-10 hour blocks of continuous good 
data in any day and to eliminate transitions into periods 
of poor atmospheric stability. On a typical day that passes 
the initial rms cut, an additional 40% of the synthesized 
data is removed due to atmospheric fiuctuations. The net 
result is that 26% of the data are kept. 

As discussed above, data are placed into 768 bins in 
right ascension with RA=0° corresponding to the first 
bin. All data falling between bins 288 and 555 are cut 
in order to eliminate observations of the Galaxy from the 
CMB analysis. This is equivalent to cutting all data with 
135° < RA< 260° centered at ^ = -60°. By observing at 
an azimuth of 208°, the galaxy cut occurs during the day, 
when the data are not generally useful for CMB observa- 
tions at 144 GHz because of the atmosphere. The Galactic 
cut overlaps well with the atmospheric cut. On a day dur- 
ing which there are minimal atmospheric fiuctuations, we 
end up cutting only 35% of the data. 

When the rms chopper position over a single chop de- 
viates by more than 0°.015 from the average position, the 
data are cut. This includes times when the chopper is in- 
tentionally shut off due to testing and maintenance as well 
as time during high winds. This cut eliminates ~ 3% of 
the data. 

The absolute timing is done through the Global Posi- 
tioning System (GPS). Data are eliminated during GPS 
drop-outs. There was also an error in the first frame of 
each logical file which made the GPS read out incorrectly 
for that frame. The first file in each set is therefore re- 
jected (« 5% of the data). The files were generally started 
directly following a calibration or during the day when we 
were at the telescope, thus the percentage of these data 
which would otherwise be retained for CMB analysis is 
small. 

11.2. Cuts to the binned data 

The goal of the cuts to the binned data is to ensure that 
only long periods of uninterrupted data are included in 
the final analysis. We cut away the low density regions by 
removing any fundamental bin with fewer than fifteen of 
the adjacent bins filled. This is done twice to ensure that 
stragglers are removed. Finally, the entire file (approxi- 
mately one day) is removed if there arc fewer than 100 
usable bins (approximately three hours of uninterrupted 
clean data). More than 90% of the days analyzed for CMB 
anisotropy have continuous 5-10 hour blocks of data. 

11.3. Offset Removal 

An offset is the value of a given harmonic when the sig- 
nal on the sky is zero. Equation 16 gives the expression 
for an element of synthesized data. Each element in the 



final data file is the average value over the RA bin, 

*n" = l^E^".i' (24) 

J = l 

where tnj is the result from each chop cycle and Ndb is 
the number of chopper cycles in a RA bin. Thus, on a 
given night, each RA bin has a single value of t*"" for each 
harmonic. These values are examined as a function of RA 
bin and both a slope and a mean are removed from each 
harmonic on each night before the nights are combined in 
order to reduce the potential effect of variations in the off- 
set over long timescales. The slope subtraction has little 
effect on the final result; its removal is prudent though not 
essential. 

Because the best data do not exactly overlap each night, 

a potential bias occurs when the mean and slope are sub- 
tracted over different regions of sky. Typical offsets for the 
TOC098 D-band data are « 150±75 in absolute value 
and drift over a timescale of days. When each full night of 
data of one harmonic is coadded and analyzed as a func- 
tion of night throughout the campaign, x^/i^ for deviation 
from a flat line is between one and four. The high-^ oflsets 
are generally more stable. Because we keep only large sec- 
tions of contiguous data and because the offset is stable, 
the effect of slightly different sky coverage per evening is 
negligible. In the analysis, we approximate the slope and 
offset removal as removing a single mean and slope from 
the entire data set as discussed below. 

The chopper-position dependent offsets (Figure 6) dis- 
cussed in Section 5 are effectively filtered by the beam 
synthesis. We check this by applying the synthesis vec- 
tors to the average chopper-synchronous signal by setting 
di = Oi in equation 16. We then examine the variation in 
the resulting synthesized offset as a function of harmonic 
and observing night. We find that the standard deviation 
of the resulting 0„ is < 20/iK for the harmonics used in the 
analysis. Because the offsets are subtracted, this variation 
does not affect the final result. 

11.4. Ergodicity of noise 

Well defined noise properties are essential for the anal- 
ysis of CMB data. After the offset subtraction and data 
cuts, we grid the i„ for each harmonic with pixel number 
vs. night of observations. We find that the average error 
bar is independent of night and pixel number. For a given 
pixel, we also check that the variation in the mean value 
from night to night is consistent with the average of the 
variances separately determined each night. 

In the high signal-to-noise channels, dT^ can be found 
by forming: 

ST^ = (AL - Al,,)/I{W) (25) 
where Aj^j is the total variance of the data for one har- 
monic and ^inst is the average variance due to instru- 
ment noise. For / < 200, the results from this calculation 
are within 10% of the results for the full likelihood anal- 
ysis (Torbet et al. 1999). At / < 150 the noise is poten- 
tially the most problematic because the offsets are gen- 
erally larger and the contribution from the atmosphere is 
larger than for the higher I data. That a simple method 
based on average noise properties gives the same answer 
(within la) as the full likelihood analysis, with its detailed 
attention to correlations and drift subtraction, gives us 



21 



Table 10 
Cuts for TOC098 D-band data 



Data cuts for the season 

Total number of days analyzed** 54 days 

Number of days used (Dl) 28 days (52%) 

Number of days used (D2) 23 days (47%) 

Percentage of data kept on good day 53% 

Net percentage kept'' 26% 



Data cuts for a typical good day Data Removed 

Galaxy 32% 

rm.s (overlaps with Galaxy cut) 24% 

Chopper 3% 

Pointing 1% 

GPS (overlaps with rms and Galaxy cut) 5% 

Net cut before beam synthesis 41% 

Net cut after beam synthesis 45% 

Net cut to the binned data'' 47% 



"This is the number of consecutive days on which the SIS system operated. After 54 days, the refrigerator cold head malfunctioned and heated 
to the point where the SISs were unusable. 

''This is the net observing efficiency for data of sufficient quality for CMB analysis. 
°This corresponds to 53% of data kept on a good day. 



confidence that the noise estimates are robust and stable. 

There is nothing in the instrument or sky of which we are 
aware that changes at the temporal frequencies and spa- 
tial scales associated with I > 200. Although the signal 
to noise is lower for large n-pt beams, the data are stable 
and insensitive to cuts. 

12. LIKELIHOOD ANALYSIS AND TESTS OF DATA 
SELECTION 

We find the most likely value of ST'^ as a function of 
angular scale with a Bayesian likelihood analysis of the 
cleaned data. The method follows the analysis of the 
Saskatoon data though the correlation matrices are con- 
siderably simpler due to the observing strategy and the 
focal plane array. One feature of the likelihood analysis 
is that that channels and harmonics can be combined in 
a precise and unambiguous way to increase the signal-to- 
noise. An added benefit of the combined bands is that 
correlation between combined bands is reduced over that 
for a single harmonic. The likelihoods for the HEMT data 
and SIS data are shown in Figure 9. 

We test that data selection eliminates the atmospheric 
and instrumental contaminants without biasing the result 
by performing a set of internal consistency tests or null 
tests. Data arc multiplied by several different sets of syn- 
thesis vectors that have no sensitivity to the sky signal and 
the analysis is carried out in the same way as for the data 
which have been multiplied by the sky-sensitive synthesis 
vectors. We determine that the data have been properly 
selected when the set of null tests fails to show a signal for 
individual synthesized beams. The data selection is there- 
fore blind to the signal on the sky, minimizing potential 
bias introduced by expectations. 

The four null tests probe different time scales. The 
quadrature signal, the diflference in signal between the first 



and second half of the chop, is sensitive to variations at 

2 Hz. For the fast and slow dither, we look at the difference 
between subsequent 0.5 second and 10 second averages re- 
spectively. We also examine the difference night to night 
and between the first and second half of the campaign. 
In the TOC098 analysis of individual harmonics for each 
channel, there was one "failure" (> 2a deviation from a 
null signal) out of 120 tests, well within expectations. The 
results of some of the null tests for the full bands are shown 
in Figure 9. 

12.1. Effects of selection criteria on the data 

Following the selection of the best cuts, we relax the se- 
lection criteria or methods in four ways. For each test, we 
perform the full likelihood analysis. We emphasize that 
the nominal cut is based on the null tests, not the data. 
When the criteria are relaxed, the null tests either show 
signs of contamination, when the likelihoods change, or 
are unaffected, when the likelihoods do not change. In 
the following we focus on the D-band data, though similar 
tests were performed on the HEMT data. 

Test 1: The minimum number of filled bins on a given 
day is increased from 100 (w 3 hours) to 250 (« 8 hours) 
and all days not common to both Dl and D2 are rejected 
from the analysis. The result is a decrease in the signifi- 
cance of the detections (due to permitting less data in the 
final analysis) but a change of less than la for the two 
low-Z points and less than a/2 for the three high I points. 

Test 2: The nominal calibration model is determined 
from fits to Jupiter and the internal calibration pulses. 
Instead, we perform the analysis simply using the calibra- 
tion value of the closest Jupiter calibration. The result 
is a change of less than a/2, primarily due to a different 
calibration. 

Test 3: We make the harmonic dependent cut to the 
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Fig. 9. — Likelihood plots for the TOCO data. All likelihoods are normalized to unity at the maximum. The x-axis is in units of 
(mK)'^. The value on the right of each panel corresponds to the effective I for each Knox filter. The thick contours are for the data. The 
thinner lines are the likelihoods for the null tests: quadrature, fast dither, slow dither, and first half minus second half signals. The / = 128 
D-band channel is possibly contaminated by the atmosphere, though the distribution of the null tests is consistent with noise. The null test 
that is close to the signal is from the slow dither and of all the null tests is the one most different from zero. Such ft: l.Scr deviations are not 
statistically unexpected. This figure clearly shows the rise and fall of a peak in the CMB angular spectrum. Furthermore, it shows that the 
null tests are consistent with zero signal. The D-band data are at I = 128, 152, 226, 306, &409. 



synthesized data more stringent. Instead of cutting 45% 
for a typical good day as shown in Table 10, 65% is cut. 
The more stringent cut is most pronounced in the lower 
harmonics because they are most affected by atmospheric 
contamination. The cut results in upper limits for the 
I = 128 and I = 158 points, consistent with the decreased 
amount of data. Above the 7-pt, the amount of data cut 
and the resulting likelihoods change very little. We also 
perform this test by eliminating just 50% (as opposed to 
45%) with similar results. 

Test 4: We apply only the initial cut to the raw data and 
so on a typical good day, 41% of the data are kept (before 
the Galaxy cut). This has the effect of introducing some 
atmospheric contamination. The result is that the I = 128 
and I = 158 points show an increase in signal of O.Su and 
1.3cr respectively while the data points corresponding to 
the highest three groups of harmonics show no significant 
change. 

We have emphasized the TOC098 D-band data because 
the points at I — 128 and / = 152 are more sensitive to 
the cuts than all the other data. This is due to the large 
atmospheric opacity at 144 GHz and large angular scales. 
The nominal cut corresponds to the best null tests and 
corresponds roughly to the sections of raw data one would 
select by eye. The D-band data for I > 200 are insensitive 
to the cuts as are all of the HEMT data. 

12.2. Correlations 

The off-diagonal terms in the noise correlation matrix, 
Cat, are small. The atmospheric correlation coefficients 
between channels (Q and Ka band) for TOC097 and be- 
tween Dl and D2 for TOC098 are of order 0.05. We ex- 
amine the autocorrelation function of the data for indi- 
vidual harmonics and find negligible correlations between 



RA bins (due to atmospheric fluctuations). The largest 
off-diagonal terms in M are < 0.4 and are in Ct- As part 
of our checks, we ensure that the likelihood is stable to 
changes of « 30 % the values of the off-diagonal terms. 

13. TESTS FOR SYSTEMATIC EFFECTS 

There are potential sources of systematic error that 
would not be revealed by the tests described above. These 
sources of error, along with their maximum contribution 
to the data, are summarized in Table 11 and described 
below. Again, we focus on the D-band data. We note that 
no effect significantly affects the final results and in par- 
ticular no effect can produce the decrease in power that is 
observed between I « 200 and I « 400. 

The systematic effects involve either beam smearing 
(Sections 13.1 - 13.5) or another mechanism (Sections 13.6 
- 13.9). For those effects related to beam smearing, the 
uncertainties are found through computing I{W) for the 
smeared beam. The fractional error for an individual har- 
monic is given by 



^ 1 d{5Ti) ^ 6./IW) 

6T, 



, , (26) 

The fractional error in groups of harmonics is smaller be- 
cause the errors are computed for the highest harmonic in 
a band group, the one with the fewest number of physical 
beams between nulls in the synthesized beam and there- 
fore the one most affected by beam smearing. 

Errors are computed for Dl, unless otherwise noted, be- 
cause it has a smaller beam than D2 and is consequently 
more affected by smearing. Nominal values of y^I{W) are 
the following: ^/I(W) = 0.377 for Dl and ^/I(W) = 0.350 
for D2 for the case of the 16-pt and 17-pt beams respec- 
tively. 
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Systematic Effect 


Maximum Resultant Change in 5T/T 


dhonner iitter 


< 0.1% 




\ U.Zi /u 


Miscalibration of chopper amplitude 


< 3% 


Mis-calibration of azimuth 


< 0.1% 


Mis-calibration of elevation 


< 1% 


Incorrect determination of beam size 


< 1% 


Electronic roll-off (between data points) 


< 0.2% 


Gain variation in offsets 


< 3% 


Chopper phase offset 


< 0.5% 



Table 11 
Systematic Effects 



13.1. Chopper Jitter 

Jitter in the chopper position, due to wind or electronic 
malfunction, smears the beam. All data for which the rms 

of the chopper position deviates more than 0.015° from the 
nominal position are rejected from the analysis. The max- 
imum resulting error in either Dl or D2 in any harmonic 
(resulting from an increase in beam size in the direction 
parallel to the chop of 0.015°) is 6^I{W) = -3.0 x lO"'^, 
or STi/Ti « 0.1%. 

13.2. Chopper Zero Offset Stability 

Beam smearing can be caused by a drift in the electronic 

zero point of the chopping mirror. The zero point is mon- 
itored throughout the campaign. We find that from the 
beginning to the end of the analyzed data the zero drifts 
by « 0°.03. We can place an upper limit on the extent to 
which this can affect our result by considering a smearing 
of the beam in the direction parallel to the chop of 0°.03. 
This results in 5s/I{W) = -7 x 10"^, or 5Ti/Ti « 0.2%. 

13.3. Mis- calibration of Azimuth 

If the azimuth changed over a period of time, the true 
beam could be uncertain by an amount equivalent to this 
change. The relative error in azimuth is < 0°.01 which 
could manifest itself as an increase in beam size in the 
direction parallel to the chop by a maximum of 0°.01. 
The resultant errors are 5^JI{W) = —2.0 x 10~^, or 
STi/Ti < 0.1%. 

13.4. Mis- calibration of Elevation 

Similarly, if the elevation changed over time, the true 
beam could be uncertain by an amount equivalent to this 

error. The relative error in elevation is < 0°.01 which could 
manifest itself as an increase in beam size in the direction 
perpendicular to the chop by a maximum of 0°.01. The re- 
sultant errors are Sy/I{W) = -3x IQ-^, or STi/Ti « 0.8%. 

13.5. Incorrect Determination of the Beam 

If the beam was incorrectly measured, there would be 
an effective beam smearing in either the direction parallel 
or perpendicular to the chop. Errors in the 6fwhm of the 
beam are « 0°.005 and « 0°.008 for Dl and D2 respec- 
tively. The maximum resulting error in either channel in 



any harmonic is 5^I{W) = 2.5 x lO'^, or 5Ti/Ti k 0.7%. 
This is an upper limit because errors in the determination 
of the beam affect the calibration as well and are therefore 
partially compensated. 

13.6. Mis- calibration of Chopper Amplitude 

If the chopper amplitude were mis-calibrated, the as- 
sumed number of physical beams fitting into a given syn- 
thesized beam could be wrong. We know the amplitude of 
the chopper to < 2% uncertainty {2a). If the value used 
in the analysis were wrong by this amount, we find the re- 
sulting change, STi/Ti < 3% for the worst case harmonic. 
Most harmonics show a change of STi/Ti < 1%. The entire 
spectrum also shifts a small amount in I. For a 2% increase 
in the chopper amplitude, we find w 2% shift in each point 
to lower I values. The sign of this effect is reversed if the 
chopper amplitude is smaller than the assumed value. 

13.7. Electronic Roll-off 

The likelihood contours shown in Figure 9 correspond 

to different post-detection frequencies as well as different 
The point at / = 226 contains data from the 6-pt to 
7-pt beams in Dl and the 6-pt to 9-pt beams in D2 cor- 
responding to / « 25 Hz. The point at I = 409 contains 
the 12-16-pt beams in Dl and 12-pt to 17-pt beams in D2, 
corresponding to / w 50 Hz. By calculation and measure- 
ment we show the electronic transfer function changes by 
< 0.2% over this range and thus does not affect the results. 

13.8. Gain Variation in the Offsets 

Variations in the cryogenic temperatures and the tem- 
perature of the warm electronics box lead to variations in 
the gain. If the effects of the decrease in gain of the warm 
amplifiers, the decrease in sensitivity of the SIS mixers, 
and the increase in output of the calibration source with 
increasing temperature are combined during the period of 
CMB observations, there could be an undetected gain drift 
of « 1%. In addition, there could be other gain effects of 
order 3% that would escape detection because the vari- 
ation in atmospheric opacity would mask their signature 
in the cal pulse. Because of the large heat capacity, such 
variations would occur on the time scale of hours. 

Changes in gain affect the offset. Since the offsets are 
different for different synthesized beams, it is possible to 
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affect the spectrum tfirough a gain change. When the av- 
erage offsets for the combined harmonics are multiphed by 
a gain variation of w 3%, a signal with amplitude w 3/iK 
at low harmonics and < 2/xK at higher harmonics results. 

Not only is this signal small, but it is largely accounted 
for through the offset and slope subtraction. 

13.9. Chopper Phase Offset 

If the phase of the chopper were to drift, the sky signal 
would be smeared because the beam would not be posi- 
tioned according to the nominal chopper template. A cut 
is made that eliminates from the analysis all data for which 
the chopper amplitude differs from the rms average value 
by more than 0.015°. This corresponds to a phase shift of 
0.5 samples. A shift of this size contributes < 0.5 /iK in 
STi. 

14. FOREGROUND EMISSION AND RESULTS 

After the raw data have been reduced and binned on 
the sky, we determine the foreground emission contribu- 
tion from our galaxy and from other galaxies. The fore- 
ground contribution at these frequencies, galactic coordi- 
nates, and angular scales is small (Tegmark et al. 2000; 
Coble et al. 1999; Masi et al. 2001). This is clear from 
Figure 10. The amplitude of the first peak is measured 
to be roughly the same from 30 to 144 GHz when the 
flux is expressed as changes in a 2.725 K thermal emit- 
ter. The spectral index of the fluctuations is (3c mb = 
ln(5Ti44/5r36.5)/lii(144/36.5) = -0.04±0.25, where ST144 
is the weighted mean of the two highest points for the D- 
band data and 6X^^,5 is a similar quantity for the HEMT 
data (36.5 GHz is the average HEMT frequency). For the 
CMB, PcMB = 0. For dust, (3b.j = 1-7 corresponds to 
PcMB = 2.05; for free-free emission Pju = —2.1 corre- 
sponds to (3c MB = —1.75. 

Fits to foreground templates were not done for Tor- 
bet et al. (1999) and Miller et al. (1999) and there 
was the possibility that the amplitude of the peak 
had a contribution from foreground emission (Page 2000; 
Knox & Page 2000). However, the mean frequency spec- 
trum of the peak clearly singles out the CMB as the 
dominant source of the fluctuations. In addition, the an- 
gular spectra of the foreground emission is much differ- 
ent than that of measured angular spectrum of the CMB 
(Miller et al. 1999). 

To quantify the foreground emission, we fit to four tem- 
plates: the SFD dust map (Schlcgel et al. 1997) {Tsfd), 
the Haslam synchrotron map (Haslam et al. 1982) (Th), 
and the radio (Tr-pt) and IRAS (Tir-pt) source com- 
pilations from the WOMBAT compilation (WOMBAT 
2001). We have not yet fit to the H^ maps that 
trace microwave free-free emission. However, we note 
fits to SK (Gaustad et al. 1996; Simonetti et al. 1996) 
found that the free-free contribution was negligible. 
In addition, cross-correlations between the WHAM Hq, 
maps (Haffncr et al. 1999) and other CMB maps, in- 
cluding QMAP, do not show a significant contribution 
(de Oliveira-Costa et al. 2001). 

There is a correlation between CMB maps at fre- 
quencies < 90 GHz and dust maps (Kogut 1996; 

This simple treatment, where we ignore chance alignments of the 
subtract the foreground in quadrature, is sufficient because the foreg 
address the full treatment and the galactic plane. 



dc Oliveira-Costa et al. 1997; Leitch et al. 1997). The 
mechanism responsible for the correlation is not 
yet certain (de Oliveira-Costa et al. 2001) though spin- 
ning dust grains is a strong candidate (Jones 1997; 
Drain & Lazarian 1998). The dust-correlated forcgroimd 
component is larger than the free-free component traced 
by H„ between 20 and 40 GHz and is not correlated to 
the Haslam map (de Oliveira-Costa et al. 1998). Thus, 
the SFD map is a good tracer of foreground emission for 
the HEMT data. The galactic latitude of the CMB scan 
covers between b = -8°, I = 280° and b = -34°, I = 335° 
as shown in Figure 8. Our sky scan passes near the Large 
Magellanic Cloud (LMC, b = -32?9, I = 280?5) and we 
remove data near it. During the day, we observe known 
sources in the galactic plane (Puchalla et al. 2001). 

The template fitting procedure is restricted to the sec- 
tion of sky analyzed for CMB anisotropy. The goal is to 
assess the contribution to the CMB results rather than 
to characterize foreground emission. First, the template 
map is multiplied by the synthesized beam in a manner 
such that the full CMB observing strategy is reproduced. 
Except in the case of the Haslam map for the I > 200 D- 
band data (where no correlation is expected), the size of 
the synthesized beam lobes is larger than the resolution of 
the template map. Next, for each harmonic, we fit to the 
template with 



= X] {tk-{aoTr-pt+aiTir-pt+a2TH+a3TsFD)) Wk 

fc=i 

(27) 

where Nn,bins is the number of filled bins around the sky 
for the n—pt function (< A^ra, tk is the CMB temperature 
(equation 24) , and Wk is the statistical weight of each mea- 
surement. The fits are performed as a function of angular 
scale and frequency. Where the fit result is significant, as 
shown in Table 12, we subtract the foreground component 
in quadrature from 5Ti as determined from the likelihood 
analysis^^. 

The fit results may be summarized as follows: 1) 
The dominant source of contamination at 30 and 40 
GHz is from the radio point sources as traced by 
the 4.85 GHz PMN catalog (GrifBth & Wright 1993; 
Condon et al. 1993; Wright et al. 1994) on which the 
WOMBAT maps are based. The extrapolation to 30 GHz 
is known to be problematic because the source spectra vary 
from source to source (Puchalla et al. 2001). We cannot 
rule out a contribution from sources that are not traced 
by the extrapolated PMN survey but, because of the mea- 
sured spectrum of the peak, such a contribution cannot be 
large. Over our sky coverage « 600 deg^, we account for 
« 100 sources. Most of the contribution comes from the 
largest 10% and is spread throughout the observing region. 
2) Contributions from synchrotron emission are negligible 
unless the synchrotron spectral index varies so much that 
the Haslam map is not a good template. The largest rms 
of the Haslam map, after applying the beam synthesis, is 
0.6 K. When extrapolated to 30 GHz with a spectral index 

CMB with the foregrounds [e.g., (de Oliveira-Costa et al. 1998)] and 
round contribution to the CMB < 4% in all cases. Future work will 



25 



Table 12 

Summary of Foreground Contributions" 



Foreground / = eO*" / = 85 / = 115 Z = 150 Z = 200 

IRAS sources at 30 GHz (^K) • • • 19 

Radio sources at 30 GHz (/xK) • • • 16 23 16 13 

Radio sources at 40 GHz (fiK) ••• 11 10 

SFD 30 GHz (pK) 18 

SFD 40 GHz liiK) ■■■ 9 



"Wc do not give any value if tlie fitted foreground contribution is < 7 /uK. At i = 60 tliis corresponds to a 1% correction and at I = 200 a 0.5% 
correction. 

''The foreground contribution for a typical value of I. 



of /3fi,7 = —2.7, one gets 5.5 fiK, consistent with the fitted 
values. The fitted dust contribution at the 30 GHz I = 60 
and 40 GHz I ~ 80 points, is small but consistent with the 
picture of 30 GHz dust-correlated emission. 3) There is no 
significant contamination at 144 GHz from either dust or 
point sources. 

Table 13 gives the results from TOGO along with the 
revised results from QMAP and SK^^ based on foreground 
analyses and new information on the calibration. The fi- 
nal results for all three experiments are shown in Figure 10 
along with a comparison to recent experiments. 

15. DISCUSSION 

Since the discovery of the anisotropy (Smoot et al. 1992), 
there have been many CMB anisotropy experi- 
ments at / < 1000 in addition to the ones 
noted thus far [ARGO (dc Bcrnardis et al. 1994; 
Ratra et al. 1999), MSAM (Wilson et al. 1999), 
UCSB/SP (Gundersen et al. 1994), White Dish 
(Ratra et al. 1998), Python (Ruhl et al. 1995), BAM 
(Tucker et al. 1996), lAB (PicciriUo & Cahsse 1993), 
Tenerife (Hancock et al. 1994), MAX (Lim et al. 1996)]. 
The SK experiment was the first to show the rise of the 
GMB angular spectrum in the region between I = 50 
and I = 250 independent of any other experiment. At 
the time, the amplitude of the peak was considered high 
though consistent with the subsequently popular A-CDM 
models. Until TOCO, the SK spectrum was in mild con- 
flict with other experiments. 

Though a straightforward read of the data prior to 
TOCO suggested a peak at I « 200 [e.g., SK and OVRO 
(Leitch et al. 1998) or SK and CAT (Baker et al. 1999)] 
there were lingering questions of cross calibration and 
point source subtraction. In fact, if one did not include 
the SK data, the combined analysis of UCSB/SP, ARGO, 
MAX, White Dish, and SuZIE favored an open universe 
(Ratra et al. 1999b). Dodelson and Knox (2000) showed 
that TOCO alone singled out a flat universe as the best 
model. Others (BahcaU et al. 1999; Bond et al. 2000), in- 
cluding Dodelson and Knox, showed that the combination 
of all the data singled out a flat universe. 

QMAP was designed to measure the anisotropy by di- 
rect mapping. Degree-resolution high signal-to-noise maps 
were made that could be compared to each other and, be- 

^® In addition to some authors on this paper, SK was analyzed by Ed 



cause they covered the same region of sky, confirmed the 
/ < 150 Saskatoon results (Xu et al. 2000). The strategy 
of direct total power mapping was later employed by the 
BOOMERanG and MAXIMA experiments, which ushered 
in very high signal-to-noise sub-degree resolution mapping. 

The TOCO experiment showed both the rise and fall 
of the first peak and showed its spectrum was that of the 
CMB with a single instrument (Figure 10). These results 
have since been amply verified by BOOMERanG, MAX- 
IMA, and DASI. The rise to the maximum has high signal- 
to-noise; the fall for / > 200 is also clearly evident though 
worth quantifying. Miller et al. (1999) reported that the 
95% upper limit of the last point in D-band, ST40Q, was 
just below the 2(t lower bound of the point of the peak, 
^2226- The penultimate point was not included in the as- 
sessment. Here we perform a more complete analysis. We 
focus on just the last three D-band points, 6T226, ST^qq, 
(5T409, because the calibration uncertainty is common to 
all and drops out of the analysis. There arc two ways a 
peak could be detected: ST226 > ST^Qg, and (5T306 > <5T409 
but with ST/^QQ < 5T22Q- The net probability that a peak 
has been detected is given by the sum of the probabilities 
of these two possibilities. 

We Monte Carlo the likelihood distributions to deter- 
mine the above probabilities and to investigate the effects 
of correlations. We find P{6T22& > ^^409) = 0.99614 and 
PiST^ofi > (5T409|(^T409 < 5T22&) = 0.00369. Thus, the 
net probability that a peak has been detected is 0.99983, 
or loosely speaking a > 3.7cr detection. The correlations 
among these three points are positive and small, of order 
0.01. When the correlations are accounted for, the net ef- 
fect is to increase the significance of a detection of a peak. 
We emphasize that the detection is model independent, 
calibration independent, and conservative in the sense that 
if there is slight contamination by point-sources or there 
is some undetected source of correlation, the probability 
of the detection increases. 

There arc now many analyses that extract cosmo- 
logical parameters assuming that the models from 
CMBFAST (Seljak & Zaldarriaga 1998) describe Na- 
ture, [e.g., (Bartlett et al. 1998; Lineweaver 1999; 
Tcgmark & Zaldarriaga 2000; i&Ec et al. 2001)]. Here 
we compute (JTpeofe and Ipeak directly from the data. This 
allows parameter estimation with the minimal amount 

WoUack and Norm Jarosik. 
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Table 13 

Summary of results from SK, QMAP, and TOCO. 



Exp eriment 


1 

te 




•» bins 






oil 




TOCO 


631^8 


2 


16,20 




35.4 


35.lt^Y 


1 9Q9+820 


QMAP 








47l? 


46.0 


48.3l?t 


240llf,« 


TOCO 




3 


17-28 


45.3t^:° 


43.4 


43.0t^:^ 


18461^^^ 


SK 




15 


24,48 


49l« 


48.0 


50.2tli 


2520l^°25 


QMAP 


1 1 1 +64 






521^ 


52.0 


54.61^:^ 


30691^5^^ 


TOCO 




6 


22-42 


70 1 +®-^ 


68.0 


67.31^:^ 


45291?^? 


QMAP 








59+? 


58.0 


60.9t?:^ 


381918?? 


TOCO 




1 


84 


54.6ti«:* 


54.6 


rq 7+18.1 
JJ- ' -16.3 


28841^72 


TOCO 


152^28 


3 


84 


82.0tii° 


82.0 


80.6llgj 


6497lli^ 


1 yJ y \J 


lOo 23 


D 


ZiJ-OO 




o i .Z 


OD.4_y ]^ 




SK 




10 


48,96 




67.6 


7o.5t^:i 


4970l^r 


TOCO 




11 


54-82 


84.7t?-i 


83.7 


82.9i?:^ 


687211JJ8 


TOCO 




6 


125 


83.0t^:« 


83.0 


81.6+?-^ 


6659l|J| 


SK 


237+11 


4 


48,96 


851^° 


83.3 


86.8li°3* 


7535_j^372 


SK 


286tlt 


4 


48,96 




84.3 


87.9tlg.5 


'7'70f;+2354 
' '^"-1720 


TOCO 




6 


165 


7o.oti;:« 


70.0 


68.8t?o'8 


4733t}^« 


SK 


349+lt 


5 


48,96 


691^^8 


67.6 


70.411^:? 


49561^3^?° 




409tg 


9 


250 


< 67 (95%con/) 




r,q q + 22.4 

ZO.O„22.4 


545 ± 2043 



"le is computed from the window function or the Knox filter following Bond (1994). In practice, we find little difference between the combined 
weighted windows (Netterfield et al. 1997) and Knox filters (Knox 1999). 

''The total number of individual n-pt functions combined in the covariance matrix M. We emphasize that all known correlations are accounted 
for. For QMAP, n-pt functions correspond to different eigenmodes for the map. 
'^The number of data points or RA bins. When two numbers or a range are given, not all n-pt beams have the same RA bins. 
''The originally published values following the convention in Netterfield et al. (1997). Calibration error is not included. 

■^The original value corrected for foreground emission (/c) . For SK, (de Oliveira-Costa et al. 1997) found a 2% contamination of the data 
due to a foreground component correlated with dust emission. For QMAP a similar correction was found (de Oliveira-Costa et al. 1999). The 
entries show the spectrum after the corrections. The corrections are done separately for Ka and Q bands before they are combined. 

^The foreground and calibration corrected values. The recalibration is based on new information about the calibration sources or, in the case 
of TOCO, on a 1-2% correction for the electronic bandpass. After SK data were published, (Mason et al. 2000) reported a more accurate 
calibration of Cas-A, which we used for QMAP. This led to an increase of 4% for the SK data and a reduction in the calibration error from 
14% to 10%. There is also a small correction to put the results in the l{l + l)/27r format as opposed to the l{2l + 1)/4it. Calibration error is 
not included. 

3 Originally, the 95% upper limit and likelihood curves were given for this bin whereas Icr error bars were plotted for the other points. We here 
adopt the convention in Mauskopf et al. (2000), and give all results with la error bars. The data are the same as before (Miller et al. 1999) 
except for the 1.7% calibration correction. The likelihood distribution is not Gaussian. The value of 23 it 22 fiK matches the distribution in 
the sense that the likelihood peak is at 23 /xK and 95% of the probability is < 67 = 23 -I- 2cr ^K. The value of 545 it 2043 (/^K)^ comes from 
fits to the likelihood following Bond et al. 2000 and is often used to represent the likelihood. Note that \/545 + 2a = 68 fiK. We set the error 
bar on STi by forming (^545 + 2ct — \/545)/2. fiK. For detailed analyses, the full likelihood as shown in Figure 9 should be used. 




Fig. 10. — The combination of the SK, QMAP, and TOCO data with all corrections. Following Netterfield et al. (1997), we plot the 
position of the maximum of the likelihood with error bars that encompass 68.3% of the likelihood for all points where there is a significant 
detection. For the last TOCO point we plot a representation of Icr based on Bond et al. (2000). The plot is made with ST as opposed to 
ST-^ because ST directly shows the signal-to-noise. For example, ST = 5 it 1 fiK has a S/N = 5 whereas the equivalent ST = 25 it 10 ^K has 
an apparent S/N ft: 3. The models are "standard CDM" (lower) and the best fit from Jaffe et al. (2000) (upper). Calibration error is not 
included. 



of model dependence. From the TOCO and SK data, the average amphtude of the peak between I = 160 and 
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/ = 240 is STpeak = 80.9 ± 3.4 ±5.1 fiK where the first 
error is statistical and second error is for calibration (the 
1 iiK shift from (Miller et al. 1999) is mostly due to the 
new calibration). Following Knox and Page (2000), we 
find Ipeak = 216 ± 14 in a relatively model independent 
parametrization 

The SK, QMAP, and TOCO experiments used a va- 
riety of techniques, separate data reduction and analysis 
pipelines, and two different calibrators. These different ex- 
periments, rich with consistency checks, trace out a peak. 
We have described the systematic checks, focusing on the 
TOCO data, and have not found any instrumental effect 
or data reduction artifact that could mimic or produce 
the signal we see. It is possible that extragalactic sources 
with spectra different from the ones we assume could alter 
our results but the effect would be small and accounting 
for it would tend to enhance the downturn for I > 220. 
In conclusion, these experiments, in particular the TOCO 
experiment, have measured the rise, amplitude, position, 
and fall of the first peak in the angular spectrum of the 
CMB. 
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